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1 Introduction 

In these lectures we describe a few applications of large iV techniques to quan- 
tum field theories (QFT) with O(N) or U(N) symmetries, where the fields are 
in the vector representation. We want to show that large iV results nicely com- 
plement results obtained from more conventional perturbative renormalization 
group (RG). Indeed the shortcoming of the latter method is that it mainly ap- 
plies to gaussian or near gaussian fixed points. This restricts space dimension 
to dimensions in which the corresponding effective QFT is renormalizable, or 
after dimensional continuation, to the neighbourhood of such dimensions. Large 
iV techniques in some cases allow a study in generic dimensions. They rely on 
noting that in the large iV limit scalar (in the group sense) composite fields have 
small fluctuations (central limit theorem). Therefore if we are able to construct 
an effective field theory for the scalars, integrating out the initial degrees of free- 
dom, we can solve the field theory in a 1/N expansion. Note that for vector 
representations the number of independent scalars is finite and independent of 
N, unlike what happens for matrix representations. This explains why vector 
models have been solved much more generally than matrix models. 

In these lectures we will in particular stress two points: first it is necessary to 
always check that the 1/N expansion is both IR finite and renormalizable. Some 
technical aspects of this question which will be described in section |4.6| . This 
is essential for the stability of the large iV results and the existence of a 1/N 
expansion. Second, the large N expansion is just a technique, with its own (often 
unknown) limitations. It should not be discussed in isolation. Instead, as we shall 
do in the following examples, it should be combined with other perturbative 
techniques and the reliability of the 1/N expansion should be inferred from the 
general consistency of all results. 

Second-order phase transitions in classical statistical physics will provide us 
with the first illustration of the usefulness of the large N expansion. Due to the 
divergence of the correlation length at the critical temperature, systems then have 
at and near T c universal properties which can be described by effective continuum 
quantum field theories. The iV-vector model that we discuss below is the simplest 
example but it has many applications since it allows to describe the critical 
properties of systems like vapour-liquid, binary mixtures, superfluid Helium or 
ferromagnetic transitions as well as the statistical properties of polymers. Before 
showing what kind of information can be provided by large iV techniques we will 
first shortly recall what can be learned from perturbative renormalization group 
(RG) methods. Long distance properties can be described in d = 4 — e dimension 
by a (4> 2 ) 2 field theory. Instead in d = 2 + e the relevant QFT model is the O(N) 
non-linear a model. It is somewhat surprising that the same statistical model 
can be described by two different theories. Since the results derived in this way 
are valid a priori only for e small, there is no overlap to test the consistency. 
The large N expansion will allow us to discuss generic dimensions and thus to 
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understand the relation between both field theories. 

Another domain of application of the large N expansion is finite size effects 
and finite temperature field theory. In these situations a dimensional crossover 
occurs between the large size or zero temperature situation where the infinite 
volume theory is relevant to a dimensionally reduced theory in the small volume 
or high temperature limit. Both effective field theories being renormalizable in 
different dimensions, perturbative RG cannot describe correctly both situations. 
Again large iV techniques will help us to understand the crossover. 

Four-fermi interactions have been proposed to generate a composite Higgs 
particle in four dimensions, as an alternative to a Yukawa-type theory, as one 
finds in the Standard Model. Again, using the specific example of the Gross- 
Neveu model, we will use large N techniques to clarify the relations between 
these two approaches. We will finally briefly indicate that other models with 
chiral properties, like massless QED or the Thirring model, can be studied by 
similar techniques. 

In the last section we return to scalar boson field theories, and examine multi- 
critical points (where the large N technique will show some obvious limitations) , 
and the double scaling limit, a toy model for discussing problems encountered in 
matrix models of 2D quantum gravity. 

2 The iV-vector model near dimension four: Renormalization Group 
(RG) 

The iV-vector model is a lattice model described in terms of iV-vector spin vari- 
ables Si of unit length on each lattice site i, interacting through a short range 
ferromagnetic O(N) symmetric two-body interaction Vij. The partition function 
of such a model can be written: 

Z= /" J] dM (S? - I) exp [-£ (S) /T] , (2.1) 

i 

in which the configuration energy £ is: 

£{S) = -Y J V* 3 Si-S J . (2.2) 

ij 

This model has a second order phase transition between a disordered phase at 
high temperature, and a low temperature ordered phase where the O(N) sym- 
metry is spontaneously broken, and the order parameter Si has a non-vanishing 
expectation value. At a second order phase transition the correlation length di- 
verges, and therefore a non-trivial long distance physics can be defined. Scaling 
and universality properties emerge which we want to study. 

To generate correlation functions one can add to £ (S) a coupling to a space- 
dependent magnetic field 

£(S) = -J2 VijSi ■ Sj + Y^Hi ■ Si. (2.3) 

ij i 
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2. 1 Mean field theory and the stability of the gaussian fixed point 

To derive the critical properties of the iV-vector model one can proceed in 
the following way: one starts from the mean field approximation, valid in high 
dimensions. One then shows that the mean field approximation is the first term 
in a systematic expansion. One discovers that for dimensions d > 4 the successive 
terms in the expansion do not modify the leading mean field behaviour. For d < 4 
instead IR divergences appear and the mean field approximation is no longer 
valid. Moreover a summation of the leading IR divergences to all orders in the 
expansion leads to an effective local (p 4 field theory. The corresponding action is 
given by the first relevant terms of the Landau-Ginzburg-Wilson hamiltonian: 



H (</>) = jd d x I c (V^) 2 + ia^(x)+6i(^ 2 (x)) : 



(2.4) 



with a, b and c being regular functions of the temperature for T close to T c . 

Note that the expression ( |2.4|) , which in the sense of classical statistical physics 
is a configuration energy, is often called hamiltonian. The reason is that if one 
starts from a classical hamiltonian and a functional integral over phase space, 
the integral over conjugate momenta is gaussian and thus trivial. From the 
point of view of quantum field theory the expression (|2.4[) has the form of an 



euclidean action, analytic continuation to imaginary time of the classical field 
theory action. We shall thus generally call it the action. 

Alternatively one can imagine starting from the configuration energy ( |2.2| ) and 
constructing Wilson's renormalization group by integrating out short distance 
degrees of freedom. The spin variable Si is then replaced by a local average, a 
vector of continuous length of the type of the field 4>(x). 

Mean field theory corresponds to the gaussian fixed point of this renormaliza- 
tion group. At the critical temperature one finds a massless free field theory 



n G {</>) = J d d xi c (v</>) : 



One then performs an analysis of the stability of the gaussian fixed point. Mean 
field theory assumes that the order parameter, here the field <fi(x), is small and 
varies only on macroscopic scales. Therefore a general action can be expanded 
in powers of the field 4>(x) and derivatives. 



J I) 



where Y^e means sum over all space integrals Tii{4>) °f O(N) symmetric monomi- 
als in <fi of degree nt and containing mi derivatives (often below called operators, 
a language borrowed from quantum field theory). 
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A convenient way to understand the relevance of the 7ii{4>) terms in the large 
distance (infrared) limit is to rescale all space or momentum variables, and mea- 
sure distances in units of the correlation length, or, at the critical temperature, 
in some arbitrary unit much larger than the lattice spacing and corresponding 
to the typical distances at which correlations are measured. 

Let us perform such a rescaling here, and rescale also the field 4>(x) in such a 
way that the coefficient of [V0(sc)j , to which all contributions will be compared, 
becomes the standard 1/2: 

x i-> Ax, (2.5) 
<P(x)^C<j)(x). (2.6) 

After this rescaling all quantities have a dimension in units of A. Our choice of 
normalization for the gradient term implies: 

C = c - 1 /2 A i-^/2 5 (2-7) 

which shows that <p now has in terms of A its canonical dimension d/2 — 1. 
A term T~Ci(<f)) then is multiplied by 

Hi{<l>) i-> A d -^ (d - 2)/2 - m ^(^>). 

For A large we observe the following: 

(i) The leading term is the term proportional to f d d x(j) 2 (x), which is multi- 
plied by A 2 . This is not surprising since it gives a mass to the field and therefore 
the theory moves away from the massless critical theory (the term is called rel- 
evant) . 

(ii) If d > 4 all other terms are multiplied by negative powers and therefore 
become negligible in the long distance limit. They are called irrelevant. The 
gaussian fixed point is stable and mean field theory thus correct. 

(iii) In four dimensions the </> 4 interaction is independent of A: it is called 
marginal while all other interactions remain irrelevant. The analysis of the sta- 
bility of the gaussian fixed point then requires a finer study which will be based 
on the field theory perturbative renormalization group. 

(iv) Below four dimensions the (f> interaction becomes relevant, the gaussian 
fixed point is certainly unstable. The question of the existence of another non- 
trivial fixed point is non-perturbative and cannot be easily answered. Partial 
answers are based upon the following assumption: the dimensions of operators 
are continuous functions of the space dimension. This means that we are go- 
ing to look for a fixed point which, when d approached four, coalesces with the 
gaussian fixed point. Moreover even at this new fixed point, at least in some 
neighbourhood of dimension four, all operators except (4> 2 ) 2 should remain ir- 
relevant. The action ( |2.4[ ) should contain all relevant operators and therefore 
enough information about the non-trivial fixed point. 
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After the rescaling (|2.5|J2.6| ) the action Ti (<p) then becomes: 



H{4>) = J d d x^[V<j ) (x)] 2 + lr<j ) 2 (x) + ^gA 4 - d ( ( j ) 2 (x)) 2 Y (2.8) 

with r = aA 2 /c, g = b/c 2 . The action ( |2.8| ) generates a perturbative expansion of 
field theory type which can be described in terms of Feynman diagrams. These 
have to be calculated with a momentum cut-off of order A, reflection of the initial 
microscopic structure. The corresponding theory is thus analogous to regularized 
quantum field theory. The precise cut-off procedure can be shown to be irrelevant 
except that it should satisfy some general regularity conditions. For example the 
propagator can be modified (as in Pauli-Villars's regularization) but the inverse 
propagator in momentum space must remain a regular function of momentum 
(the forces are short range). 

Let us call r c the value of the parameter r which corresponds, at g fixed, to 
the critical temperature T c at which the correlation length £ diverges. In terms 
of the scale A the critical domain is then defined by: 

physical mass = £ _1 C A =^> \r — r c \ <C A 2 

distances ^> 1/A or momenta <C A, (2.9) 

magnetization M = \ (<f>(x)) | < C _1 ~ A (d/2)_1 . 

Note that these conditions are met if A is identified with the cut-off of a usual field 
theory. However an inspection of the action ( |2.8|) also shows that, in contrast with 
conventional quantum field theory, the </> 4 coupling constant has a dependence 
in A given a priori. For d < 4 the ^> 4 coupling is very large in terms of the 
scale relevant for the critical domain. In the usual formulation of quantum field 
theory instead the bare coupling constant also is an adjustable parameter. This 
implies for instance that for d < 4 (super-renormalizable theory) the coupling 
constant varies when the correlation length changes. This is a somewhat artificial 
situation if one believes that that the initial bare or microscopic theory has a 
physical meaning. 

The critical properties of the field theory (like the long distance behaviour 
of correlation functions) can then be analyzed by RG methods in 4 — e dimen- 
sion, i.e. near the so-called upper-critical dimension (and with some additional 
assumptions in d < 4). 

Dimensions of fields. Because we deal with translation invariant theories, we 
will generally discuss the scaling behaviour of correlation functions in momentum 
variables. Let us relate the scaling behaviour of connected correlation functions 
expressed in terms of space and momentum variables. When functions have a 
scaling behaviour, one defines 



i[Oi(xi/\)) = X D (YlOiixi)) with D = J2d 0i , (2.10) 
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where Oi, sometimes called operator, is a local polynomial in the basic fields 
(associated with the order parameter), and the quantity do i , which we sometimes 
also denote [Oi], is called the dimension of the field (operator) Oi. 

After Fourier transformation and factorization of the 5-function of momentum 
conservation, one then finds, in d space dimension, 

(j[Oi(\ Pi \ =A D 7[[ 'W) with now D' = d + J2(do i ~d). (2.11) 

\i=l I c \ i I c i 

Finally it is convenient to introduce the Legendre transform T(<p) of the gener- 
ating functional W(H) = TlnZ of </>-field connected correlation functions. We 
denote by and the corresponding connected and 1PI functions. One 

verifies that if one performs a Legendre transformation on the source associated 
with the field (operator) Oi, the quantity do i — d in equation ( |2.11| ) is replaced 
by -d 0i - 



2.2 RG equations for the critical (massless) theory 

The field theory with the action (|2.8| ) can now be studied by field theoretical 
methods. From simple power counting arguments one concludes that the critical 
(or massless) theory does not exist in perturbation theory for any dimension 
smaller than 4. If we define, by dimensional continuation, a critical theory in 
d = 4 — e dimensions, even for arbitrarily small e there always exists an order in 
perturbation (~ 2 /e) at which IR (infrared) divergences appear. Therefore the 
idea, originally due to Wilson and Fisher, is to perform a double series expansion 
in powers of the coupling constant g and e. Order by order in this expansion, 
the critical behaviour differs from the mean field behaviour only by powers of 
logarithm, and we can construct a perturbative critical theory by adjusting r to 
its critical value r c (T = T c ). 

To study the large cut-off limit we then use methods developed for the con- 
struction of the renormalized massless (p 4 field theory. We introduce rescaled 
(renormalized) correlation functions, defined by renormalization conditions at a 
new scale \i <C A, and functions of a renormalized coupling constant g r . We write 
here equations for Ising-like systems, the field <p having only one component. The 
generalization to the N- vector model with O(N) symmetry, is straightforward 
except in the low temperature phase or in a symmetry breaking field, a situation 
which will be examined in section |2l]. Then: 

If) ( P ;g T ,iJ,A)\ p2=0 = 0, 
< ^r^(p;^,A)| p2= ^ = l, (2.12) 
(pi = fj,6i;g Y ,/j,A) = /j £ g T , 
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in which 9i is a numerical vector. These correlation functions are related to the 
original ones by the equations: 



r! B) (ft; 9r, A*) A) = Z n / 2 (g, \/ y) r< B > ( Pi - g, A) . 



(2.13) 



Renormalization theory tells us that the functions If (Pi', g T , At, A) of equation 
(|2.13|) have at pi, g r and /i fixed, a large cut-off limit which are the renormalized 
correlation functions 1; (pi;g T1 fx). A detailed analysis actually shows that at 
any finite order in perturbation theory: 



rf ) ( Pi ;g r ,fi,A) = lf»> {pi;g T ,n) + (A" 2 (lnA) L ) 



(2.14) 



in which the power L of In A increases with the order in g and e. Furthermore 
the renormalized functions Ir do not depend on the specific cut-off procedure 
and, given the normalization conditions ( |2.12p , are therefore universal. Since 

the renormalized functions Ir and the initial ones l?( n ) are asymptotically pro- 
portional, both functions have the same small momentum or large distance be- 
haviour. To determine the universal critical behaviour it is thus sufficient to 
study the renormalized field theory. And indeed most peturbative calculations 
of universal quantities have been performed in this framework. However, it is 
interesting to determine not only the asymptotic critical behaviour, but also the 
corrections to the asymptotic theory. Furthermore, renormalized quantities are 
not directly obtained in non-perturbative calculations. For these reasons it is 
also useful to express the implications of equation (|2.13j ) directly on the initial 
theory. 

Bare RG equations. Let us differentiate equation (|2.13j ) with respect to A at 
g r and \i fixed, taking into account (|2.14j ): 



A 8 

<9A 



Z n ' 2 (g, A//u) rW ( Pi ; g, A) = O (A _2 (ln A) L ) . (2.15) 



g r ,/i fixed 



We now neglect corrections subleading (in perturbation theory) by powers of A. 
Then, using chain rule, we can rewrite equation ( |2.15|) as: 



aJ^( 9 ,a/„)| 



II 



(p i ;p,A) = 



(2.16) 



The functions (3 and 77, which are dimensionless and may thus depend only on 
the dimensionless quantities g and A///, are defined by: 



V (9, A///) 



A 



d_ 



-A 



d_ 

dX 



hiZ(g,A/n) 



(2.17a) 
(2.176) 



9r,M 
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However, the functions (3 and r\ can also be directly calculated from equation 
(|2.16|) in terms of functions which do not depend on /i. Therefore the 



functions (3 and r\ cannot depend on the ratio A/ ' \x (in the definitions ( j2.17| ) 



consistency requires that contributions which goes to zero like some power of 
/(//A, should be neglected, as in equation ( 2.15| )). Then equation can be 

rewritten: 

On 

A aA +/3(9) aJ-2" (9) 



r'")(p,; 9 ,A) = 0. (2.18) 



Equation ( |2.18| ) is an equation satisfied when the cut-off is large by the physical 



correlation functions of statistical mechanics which are also the bare correlation 
functions of quantum field theory. It expresses the existence of a renormalized 
theory. 

2.3 RG equations and large distance behaviour: the e-expansion 

Equation ( |2.18|) can be solved by the method of characteristics: one introduces 
a dilatation parameter A, together with a running coupling constant g(X) and a 
scale dependent field renormalization Z(X) satisfying the flow equations 

X-^g (X) = (3 (g (X)) , g(l)=g; (2.19a) 
A-^ In Z (A) = V (g (A)), Z(l) = l. (2.196) 

The behaviour of correlation functions for | A (A — > 0) is then governed by 
IR fixed points, zeros of the RG /3-function with a positive slope. 

The RG functions f3 and n can be calculated in perturbation theory. From 
the relation between bare and renormalized coupling constant and the definition 
(|2.17a|) it follows that (e = 4 - d): 



N + 8 



(3 (g, e) = -eg + -r^g 2 + O {g 3 , g 2 e) . (2.20) 

Let us now assume that g initially is sufficiently small, so that perturbation 
theory is applicable. We see that above or at four dimensions, i.e. e < 0, the 
function (3 is positive and g(X) decreases approaching the origin g = 0. We 
recover that the gaussian fixed point is IR stable for d > 4, and find that it ia 
also stable at d = 4. 

Below four dimensions, instead, the gaussian fixed point g = is IR repulsive. 
However, expression ( |2.20| ) shows that, for e small, /3(g) now has a non-trivial 
zero g*: 

W) = 0, 9* = §^ l e + 0(e 2 ), with /3'(g*) = u = e + O (e 2 ) . (2.21) 
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The slope u at the zero is positive. This non-gaussian fixed point thus is IR 
stable, at least in the sense of an e-expansion. In four dimensions it merges with 
the gaussian fixed point and the eigenvalue u vanishes, indicating the appearance 
of a marginal operator. 

The solution of the RG equation then determines the behaviour of r( n ) (p$; g, A) 
for \pi\ <C A: 

rW{\ Pi ;g,A) ~ X d -^ 2 ^ d - 2+ ^T^( Pi ;g*,A), (2.22) 

A — >0 

where rj = Tj(g*). Critical correlation functions have a power law behaviour for 
small momenta, independent of the initial value of the </> 4 coupling constant g. 

In particular the small momentum behaviour of the inverse two-point function 
is obtained for n = 2. For the two-point function (p) this yields: 



w {2) ( P ) = \r( 2) ( P )} 1 , ~ o i /p 2 ~ v - (2-23) 



The spectral representation of the two-point function implies r\ > 0. A short 
calculation yields: 

V = ^e 2 + 0(e*). (2.24) 

The scaling in equation fl2.22p indicates that the field <f>(x), which had at the gaus- 
sian fixed point a canonical dimension (d — 2)/2, has now acquired an "anoma- 
lous" dimension (see the discussion of the end of section |2.1|) : 

= \ {d - 2 + rj) . 



These results call for a few comments. Within the framework of the ^-expansion, 
one thus proves that all correlation functions have, for d < 4, a long distance 
behaviour different from the one predicted by mean field theory. In addition 
the critical behaviour does not depend on the initial value of the tfi 4 coupling 
constant g. At least for e small one may hope that the analysis of leading IR 
singularities remains valid and thus it does not depend on any other coupling 
either. Therefore the critical behaviour is universal, although less universal than 
in mean field theory, in the sense that it depends only on some small number of 
qualitative features of the system under consideration. 

2-4 Critical correlation functions with <p 2 (x) insertions 

RG equations for critical correlation functions with J d d x(j) 2 (x) insertions can 
also be derived. The operator <p 2 (x) has a direct physical interpretation. It is 
the most singular part (i.e. the most relevant) of the energy density (|2.4j) . Long 
distance scaling properties follow. Moreover these RG equations can be used to 
derive RG equations for correlation functions in the whole critical domain. 
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We denote by r('' n ) (qi, . . . , qi; pi, . . . , p n ; g , A) the mixed 1PI correlation func- 
tions of the order parameter 4>(x) and the energy density ^(f) 2 (x) (n<p fields and 
1 7}(f) 2 operators, with (I + n) > 2). Renormalization theory tells us that we can 

define renormalized correlation functions r^'" - '' {qi',Pj; g T , fi) which, in addition to 
conditions ( |2.12| ), satisfy: 



9 2 = |M 2 



1 





(2.25) 



r^' 2) (q;pi,p 2 ;g I ,fi) 

r?' 0) (9, -?;&,//) 

and are related to the original ones by: 



lim Z n / 2 (Z 2 /Z) 1 [r^ (quPj; g, A) - 5 n0 5 l2 A~ e A\ = if^ {q i]Pj] g r , //) . 

A^oo L 

(2.26) 

Z 2 (g, A///) and A (<?, A///) are two new renormalization constants. 

Differentiating with respect to A at g T and \i fixed, as has been done in section 



2T2| , and using chain rule one obtains a set of RG equations: 

r ( '< n) = S n0 8 l2 A- £ B(g) 



On 
A 0A +(3{9) d-g-2 r]{9) -^ {9) 



(2.27) 



In addition to (3 and r\ (equations ( |2.17|) ) two new RG functions, rj 2 (g) and B(g). 
appear: 



m (g) 
Big) 



-A 



A 



d_ 

dX 

d_ 



ln[Z 2 (g,A/ri/Z(g,A/(j,)] 



2^72 (g) - e 



A(g,A/jj) 



(2.28) 
(2.29) 



Note that for n = 0, I = 2, the RG equation ( 2.27 ) is not homogeneous. This is a 
consequence of the non-multiplicative character of renormalization in this case. 

In the homogeneous case, equation ( |2.27| ) can be solved exactly in the same 
way as equation (|2.16|) . A new function ( 2 (\) has to be introduced, associated 
with the RG function rj 2 (g)- Again the solution of equation (|2.27|) combined with 
simple dimensional analysis leads to the scaling behaviour 



(\ qi ;\ Pj ;g,A) oc X d-n(d-2+ v )/2-l/u ^ 
X — >0 

where the correlation length exponent v is related to rj 2 (g*) by: 

v=[r ]2 (g*) + 2}- 1 . 



(2.30) 



(2.31) 
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The dimension of the field <fi 2 follows (see section |2.1|) 



d^ = d-l/u. (2.32) 

Using equations (|2.28| , |2.3l|) it is easy to calculate 772 (#) at one-loop order. At the 
fixed point g = g* (equation ( j2.21| )) one then obtains the exponent v: 



2, = l + gt| £ + 0( £ 2 ) 



The (^ 2 ^) 2 ) correlation function. The </> 2 (energy density) two-point function 
p(2,o) sa ti s f] e g an inhomogeneous RG equation. To solve it one first looks for a 
particular solution, which can be chosen of the form A~ e C2(g): 

d (9) C 2 (g) -[e+ 27? 2 (g)] C 2 (g) = B (g) . (2.33) 

The solution is uniquely determined by imposing its regularity at g = g*. 

The general solution of equation (|2.27 ) is then the sum of this particular 



solution and of the general solution of the homogeneous equation which has a 
behaviour given by equation ( |2.30|) : 



T^(Xq;g,A)-A- s C 2 (g) ~ \ d ~ 2 ^ . (2.34) 

A — >0 



Remarks. 

(i) The physics we intend to describe corresponds to integer values of e, 
e = 1,2. Although we can only prove the validity of all RG results within 
the framework of the ^-expansion, we shall eventually assume that their valid- 
ity extends beyond an infinitesimal neighbourhood of dimension 4. The large 
iV-expansion provides a test of the plausibility of this assumption. The decisive 
test comes, of course, from the comparison with experimental or numerical data. 

(ii) In four dimensions the </> 4 interaction is marginally irrelevant; the renor- 
malized coupling constant of the </> 4 field theory goes to zero only logarithmically 
when the cut-off becomes infinite. This induces logarithmic corrections to mean 
field theory. Moreover, since no other fixed point seems to exist, this leads to 
the so-called triviality property (see section |2.7| ) of the </> 4 quantum field theory. 

2.5 Scaling behaviour in the critical domain 

We have described the scaling behaviour of correlation functions at criticality, 
T = T c . We now consider the critical domain which is defined by the property 
that the correlation length is large with respect to the microscopic scale, but 
finite. 

Remark. The temperature is coupled to the total configuration energy. There- 
fore a variation of the temperature generates a variation of all terms contributing 
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to the effective action. However the most relevant contribution (the most IR sin- 
gular) corresponds to the <f> 2 {x) operator. We can therefore take the difference 
t = r — r c (x T — T c between the coefficient of <p 2 in ( |2.8| ) and its critical value 
as a linear measure of the deviation from the critical temperature. Dimensional 
analysis then yields the relation 

( Pi ; t, g, A) = Ad-n(d-2)/2 T (n) ( p . A -i. f A -2 ? ^ ^ _ ( 2 .35) 

With this parametrization the critical domain corresponds to \t\ C A 2 . 

Expansion around the critical theory. One thus adds to the critical action a 
term of the form wt J d d xcj) 2 (x). To derive RG equations in the critical domain 
one expands correlation functions in formal power series of t. The coefficients 
are critical correlation functions involving 4> 2 (x), for which RG equations have 
been derived in section |2.4| , inserted at zero momentum. Some care has to be 
taken to avoid obvious IR problems. Summing the expansion, one obtains RG 
equations valid for T ^ T c , \T — T c \ <C 1. 

After summation of the t-expansion one finally obtains the RG equation: 



n 



r<»>(Pi;<,S.A) = 0. (2.36) 



Scaling laws above T c . As for previous RG equations, equation ( |2.36|) can be 



integrated by using the method of characteristics. In addition to the functions 
g(X) and Z(X), one needs a running temperature t(X). Taking the large A, or 
the small A limit one finally obtains: 

r» (pi;t,g,A = l) ~ m^-^-^^F^ {pi/m), (2.37) 
lpil<l 

with: 

m(A = l) (2-38) 
From equation ( |2.37| ) we infer that the quantity m is proportional to the physical 



mass or inverse correlation length. Equation ( 2.38 ) then shows that the diver- 
gence of the correlation length £ = m _1 at T c is characterized by the exponent 
v. 

For t 7^ 0, the correlation functions are finite at zero momentum and behave 
as: 

(0; t, g, A) oc t^-^-2+r/)/2)_ (2.39) 

In particular for n = 2 we obtain the inverse magnetic susceptibility: 

x - 1 = rW(p = 0;t,g,A)Kt^ 2 - r i\ (2.40) 

The exponent which characterizes the divergence of x is usually called 7. The 
equation (|2.39| ) establishes the relation between exponents: 

7 = 1/(2-77). (2.41) 
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2.6 Scaling laws in a magnetic field and below T c 

In order to pass continuously from the disordered (T > T c ) to the ordered 
phase (T < T c ), avoiding the critical singularities at T c , it is necessary to add 
to the action an interaction which explicitly breaks its symmetry. One thus add 
a small magnetic field to the spin interactions. One then derives RG equations 
in a field, or at fixed magnetization. In this way correlation functions above 
and below T c can be continuously connected, and scaling laws established in the 
whole critical domain. The first example is provided by the relation between 
field and magnetization, i.e. the equation of state. 

The equation of state. Let us call M the expectation value of 4>(x) in a constant 
field H (for N > 1 the quantities M and H should be regarded as the length of 
the corresponding vectors) . The thermodynamic potential per unit volume, as a 
function of M, is by definition: 



Q^T (M, t, g,\) = J2 -j-rW ( Pi = 0; t, g, A) . (2.42) 
The magnetic field H is given by: 



dM n\ 

n=l 

Noting that n = M (d/dM), we immediately derive from the RG equation (|2.36|) , 
the RG equation satisfied by H (M, t, g, A): 

A M + 13 « | " & (9) (l + U±) - m to) «|] H (Af , t, g , A) = . 

(2.44) 

To integrate equation ( |2.44j ) by the method of characteristics we have to intro- 
duce, in addition to the functions (7(A), t(X) and Z(X), a new function M(X). 
However one verifies that M(X) is given by M(A) = MZ-V2(A). 

Then from the arguments outlined in previous sections one derives the scaling 
form 

H (M, t, g, 1) ~ M s f (tM-W") , (2.45) 

with: 

(3=\v{d-2 + r 1 ) = vd 4> , S= d ] + 2 ~ V = (2.46) 

Equation ( |2.45| ) exhibits the scaling properties of the equation of state. Moreover 
equations (|2.46| ) relate the traditional critical exponents which characterize the 
vanishing of the spontaneous magnetization and the singular relation between 
magnetic field and magnetization at T c respectively to the exponents n and v 
introduced previously. 
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The universal function f(x) is infinitely differentiable at x = 0. because when 
M is different from zero the theory remains massive even at t = 0. The magnetic 
field H has a regular expansion in odd powers of M for t > 0. This implies 
that when the variable x becomes large and positive, f(x) has the expansion 
(Griffith's analyticity): 

oo 
p=0 

The appearance of a spontaneous magnetization, below T c , implies that the 
function f(x) has a negative zero Xq. Then equation (|2.45|) leads to the relation: 

M=\x \~ P (-tf for H = 0, t<0. (2.48) 



Equation ( 2.48|) gives the behaviour of the spontaneous magnetization when the 



temperature approaches the critical temperature from below. 

Correlation functions in a field. We now examine the behaviour of correlation 
functions in a field. We write expressions for Ising-like systems. In the ordered 
phase some qualitative differences appear between systems which have a discrete 
and a continuous symmetry. We illustrate these differences with an example at 
the end of the section. 

The correlation functions at fixed magnetization M are obtained by expand- 
ing the generating functional V (M (x)) of 1PI correlation functions, around 
M (x) = M. From the RG equations satisfied by the correlation functions in 
zero magnetization (equations ( |2.36| )) it is then easy to derive: 



d d 1 / d \ d 

^ + (3(9)^- «V (9)(n + M— - m (g) t- 



eA—ag 2— V dMj -^^j^W^AHO. 

(2.49) 

This equation can be solved by exactly the same method as equation ( |2.44j) . One 
finds 

r( n) (p i; t, M,g, A = 1) ~ m ld-(d-2+v)m F (n) (pi/m, tm' 11 ^ , (2.50) 

for \pi\ < 1, |t| <C 1, M 1 and with the definition: 

m = M v ' p . (2.51) 

The r.h.s. of equation ( |2.50| ) now depends on two different mass scales: m = 
M v /P and t v . 

Correlation functions below T c . We have argued above that correlation func- 
tions are regular functions of t for small t, provided M does not vanish. It is 
therefore possible to cross the critical point and to then take the zero external 
magnetic field limit. In the limit M becomes the spontaneous magnetization 



17 



which is given, as a function of t, by equation ( 2.48 ). After elimination of M in 
favour of t in equation ( |2.50|) , one finds the critical behaviour below T c : 



rW ( Pi ;t, M (t, H = 0) , g, 1) ~ m d-n{d-2+ v )/2 F {n) ^ ^ m j ^ ^ 2 _ 52 j 

with: 

m = laol"" (-*)" , i? = 0, £<0. (2.53) 

We conclude that the correlation functions have exactly the same scaling be- 
haviour above and below T c . 

The extension of these considerations to the functions with 4> 2 insertions, T^ l,n ^ 
is straightforward. In particular the same method yields the behaviour of the 
specific heat below T c : 

r^(q = 0,M(H = 0,t))-A- £ C 2 (g) ~ A~ (-t)~ a , (2.54) 

for t < 

which similarly proves that the exponent above and below T c are the same. 

Note that the constant term A~ e C2 (g) which depends explicitly on g is the 
same above and below T c , in contrast with the coefficient of the singular part. 

The derivation of the equality of exponents above and below T c , relies on the 
existence of a path which avoids the critical point, along which the correlation 
functions are regular, and the RG equations everywhere satisfied. 

Remark. The universal functions characterizing the behaviour of correlation 
functions in the critical domain still depend on the normalization of physical 
parameters t, H, M, distances or momenta. Quantities which are independent of 
these normalizations are universal pure numbers. Simple examples are provided 
by the ratios of the amplitudes of the singularities above and below T c like A + j A~ 
for the specific heat. 

The O(N)- symmetric N -vector model. We now indicate a few specific proper- 
ties of models in which the action has a continuous O(N) symmetry. 

The differences concern correlation functions in a field or below T c . The addi- 
tion of a magnetic field term in an O(N) symmetric action has various effects. 

First, the magnetization and the magnetic field are now vectors. The RG 
equations have exactly the same form as the Ising-like N = 1 case but the 
scaling forms derived previously apply to the modulus of these vectors. 

Second, since magnetic field or magnetization distinguish a direction in vector 
space, there now exist 2 n n-point functions, each spin being either along the 
magnetization or orthogonal to it. When the continuous 0(N) symmetry of the 
action is broken linearly in the dynamical variables (as in the case of a magnetic 
field) these different correlation functions are related by a set of identities, called 
WT identities. The simplest one involves the 2-point function T^ 1 , at zero 



18 



momentum, of the components orthogonal to M, i.e. the transverse susceptibility 

T t ( P = 0)=Xt 1 =H/M. (2.55) 

It follows that if H goes to zero below T c , H/M and therefore Tt at zero mo- 
mentum vanish. The latter property implies the existence of N — 1 (massless) 
Goldstone modes corresponding to the spontaneous breaking of the O(N) sym- 
metry. 

Note finally that the inverse longitudinal 2-point function T-^(p) has IR sin- 
gularities at zero momentum in zero field below T c generated by the Goldstone 
modes. This is characteristic of continuous symmetries, and will play an essential 
role in next section. 

2.7 Four dimensions: logarithmic corrections and triviality 

Let us just briefly comment about the situation in four dimensions. If we solve 
the RG equation 

A-^<7(A) = /?(<7(A)), 

for the running coupling constant, assuming that (3(g) remains positive for all 
g > (no non-trivial fixed point), we find that g(X) goes to zero logarithmically; 
the operator </> 4 is marginally irrelevant. Writing generally 

P(g) = P29 2 + p39 3 + 0(g 4 ), /3 2 >0, 

we find for A — > 0: 

lnA= -/d<A)-| ln9(A)+x(9) ' (2 - 66) 

with: 

A ' (9) = h ~ I ln9 ~ f d9 ' (m + w) ■ 

Since the running coupling constant goes to zero in the long distance limit, 
quantities can be calculated from perturbation theory. From the point of view of 
critical phenomena logarithmic corrections to mean field theory are generated. 

Finally let us note that empirical evidence coming from lattice calculations 
strongly suggests the absence of any other fixed point. 

From the point of view of particle physics one faces the triviality problem: for 
any initial bare coupling constant g the renormalized coupling g(/j,/A) at scale \x 
much smaller than the cut-off A behaves like 
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(2.57) 



Therefore if one insists in sending the cut-off to infinity one finds a free (trivial) 
field theory. However in the modern point of view of effective field theories, one 
accepts the idea that quantum field theories may not be consistent on all scales 
but only in a limited range. Then the larger is the range the smaller is the low 
energy effective coupling constant. In the standard model these comments may 
apply to the weak-electromagnetic sector which contains a (j) 4 interaction and 
trivial QED. 
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3 The O(N) Spin Model at Low Temperature: the Non-Linear cr-Model 



Let us again consider the lattice model (|2.1| , |2^ ) of section ^ with partition func- 
tion: 

Z = J J] dS, 5 (S 2 - 1) exp J2 V *J S * ■ Sj/ t ■ 

i ij 

We will now discuss this model from the point of view of a low temperature 
expansion. The methods we employ, however, apply only to continuous symme- 
tries, here to N > 2. They rely on the property that models with continuous 
symmetries, in contrast to models with discrete symmetries, have a non-trivial 
long distance physics at any temperature below T c , due to the massless Goldstone 
modes. 

We first prove universal properties of the low temperature, ordered, phase 
at fixed temperature. Then, in the non-abelian case, iV > 2, we show that 
additional information about critical properties can be obtained, by analyzing 
the instability of the ordered phase at low temperature and near two dimensions, 
due to Goldstone mode interactions. 

The analysis is based on the following observation: The iV-vector model 
( |2.1| , |2.2| ) can be considered as a lattice regularization of the non-linear a-model 
(note 2Si • S 3 • = 2 — (Sj — Sj) 2 ). The low temperature expansion of the lattice 
model is the perturbative expansion of the regularized field theory. The field 
theory is renormalizable in dimension two. RG equations, valid in two and more 
generally 2 + e dimension follow. Their solutions will help us to understand the 
long distance behaviour of correlation functions. 

It is somewhat surprising that two different continuum field theories, the (4> 2 ) 2 
and the non-linear cr-model describe the long distance physics of the same lattice 
model. This point will be clarified by an analysis of the 1/iV-expansion of both 
field theories. This property, totally mysterious at the classical level, emphasizes 
the essential nature of quantum (or statistical) fluctuations. 

3.1 The non-linear a-model 

We now study the non-linear c-model from the point of view of renormalization 
and renormalization group. In continuum notation the field S(x) has unit length 
and the action is 



where t is proportional to the temperature T. To generate perturbation theory 
we parametrize the field S(x): 

S(x) = {cr(x),7r(x)}, 

and eliminate locally the field cr(x) by: 

a(x) = (1-tt 2 (x)) 1/2 . 
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This parametrization is singular but this does not show up in perturbation theory 
which assumes tt(x) small. 

The O(N) symmetry. The 0(N — 1) subgroup which leaves the component a 
invariant acts linearly on the N — 1 component vector n. However a general O(N) 
transformation will transform it into a linear combination of n and a/1 — tt 2 - 
The O(N) symmetry is realized non-linear ly. An infinitesimal transformation 
corresponding to the generators of O(N) not belonging to 0(N — 1) takes the 
form 

5n = ojs/l — 7T 2 , 

where u is a N — 1 component vector of parameters corresponding to these 
generators. 

As we have done for the (</> 2 ) 2 model, we scale all distances in order to measure 
momenta in units of the inverse lattice spacing A. We thus write the partition 
function: 

Z = J [(l-7r 2 (x))" 1/2 d7r(x)]exp[-5(7r)], (3.1) 



with 



A d ' 2 f 

S ^ = l ddx 9ij( n ) d ^i( x ) d ^j( x )^ (3-2) 

where gij is the metric on the sphere 

9ij = Sij + J ri7r \ - (3-3) 

1 — 7T^ 

Moreover, as expected, the functional measure is related to the metric by 

\/ det ^ = 7f=f • 

Propagator, perturbation theory and power counting. Unlike the (j) 4 field theory, 
the action is non-polynomial in the fields. An expansion of the action in powers 
of 7r generates an infinite number of interactions. However we note that the 
power of t in front of a diagram counts the number of loops. Therefore at a finite 
loop order, only a finite number of interactions contribute. 

The 7T propagator is proportional to: 

A w (fc) = 

The 7r thus has the usual canonical dimension (d — 2)/2. Since we have interac- 
tions with arbitrary powers of tt, the model is renormalizable in two dimensions, 
where all interactions have dimension two. 
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The role of the functional measure. If we try to write the functional measure 
as an additional interaction we find 

II h 1 m , =expHEMi-^))- 
x v 1 -*^) x 

This quantity is well-defined on the lattice but not in the continuum. This prob- 
lem, which already appears in quantum mechanics (d = 1) reflects the necessity 
of a lattice regularization to precisely define the quantum hamiltonian in the 
presence of interactions with derivatives. A perturbative solution is provided by 
dimensional regularization, where this term can simply be omitted. In lattice 
regularization it cancels quadratic divergences. 

IR divergences, spontaneous symmetry breaking and the role of dimension two. 
We see that the perturbative phase of the non-linear a model is automatically 
a phase in which the O(N) symmetry is spontaneously broken, and the (N — 1) 
components of S(x), ir(x), are massless Goldstone modes. 

(i) For d < 2 we know from the Mermin-Wagner theorem that SSB with order- 
ing ((S) ^ 0) is impossible in a model with continuous symmetry and short range 
forces. Correspondingly IR divergences appear in the perturbative expansion of 
the non-linear a model for d < 2, for example (a) diverges at order t as f d d p/p 2 . 
For d < 2 the critical temperature T c vanishes and perturbation theory makes 
sense only in presence of an IR cut-off which breaks explicitly the symmetry and 
orders the spins (thus selecting a classical minimum of the action). Therefore 
nothing can be said about the long distance properties of the unbroken theory 
directly from perturbation theory. 

(ii) For d > 2 instead, perturbation theory which predicts spontaneous sym- 
metry breaking (SSB), is not IR divergent. This is consistent with the property 
that in the iV-vector model, for d > 2, the O(N) symmetry is spontaneously 
broken at low temperature. At T < T c fixed, the large distance behaviour of the 
theory is dominated by the massless or spin wave excitations. On the other hand 
nothing can be said, in perturbation theory, of a possible critical region T ~ T c . 

To go somewhat beyond perturbation theory we shall use field theory RG 
methods. It is therefore necessary to first define the model in two dimensions 
where it is renormalizable. There IR divergences have to be dealt with. We 
thus introduce an IR cut-off in the form of a magnetic field in the a direction (a 
constant source for the a field) 

-hy/1 -7r 2 (:r)J • 

(3.4) 

Expanding the additional term in powers of % we see that it generates a mass 
term 



S(7T,h) 



(tT ■ dfjTtjx))' 
1-7V 2 (X) 
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and additional interactions of dimension in d = 2. 

We then proceed in formal analogy with the case of the (cj) 2 ) field theory, 
i.e. study the theory in 2 + e dimension as a double series expansion in the 
temperature t and s. In this way the perturbative expansion is renormalizable 
and RG equations follow. 

3.2 RG equations 

Using power counting and some non-trivial WT identies (quadratic in the 1PI 
functional) one can show that the renormalized action takes the form: 



<S r (7T r , h T ) 



d-2 ■ 



2t r Z t 



d d x 



(<9 M 7r r ) 2 + (d^a r ) 2 



d-2 



in which [i is the renormalization scale and: 



K \ ar(x)d d x, (3.5) 



(3.6) 



Note that the renormalization constants can and thus will be chosen h indepen- 
dent. This is automatically realized in the minimal subtraction scheme. 
The relation: 

Tr r (x) = Z-^ 2 tt(x), (3.7) 



implies 



fJL d-2h L = A d-2 z l/2^ (gg) 

t 

With our conventions the coupling constant, which is proportional to the tem- 
perature, is dimensionless. The relation between the cut-off dependent and the 
renormalized correlation functions is: 



Z n ' 2 (A/p, t) rW ( Pi ; t, h, A) = rW ( Pi ; t r , h r , 



(3.9) 



Differentiating with respect to A at renormalized parameters fixed, we obtain 
the RG equations: 



r (B) fat, h,A) = o, 



(3.10) 



where the RG functions are defined by: 



Pit) = 

««> = A Jv 



t, 



ren. fixed 



(-InZ) 



(3.11) 



ren. fixed 



ln/i. 



ren. fixed 
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The coefficient of d/dh can be derived from equation ( |3.8|) which impiies (taking 
the iogarithm of both members): 



= h-'A-^h + d-2- - ^f-, (3-12) 

and therefore: 

p(t) = 2-d+ 1 -((t) + ^. (3.13) 

To be able to discuss correlation functions involving the a-field, we also need the 
RG equations satisfied by connected correlation functions W^: 



. d „, , d n .. s /L, , Bit) \ , d 
A 7TT + W-ST + n((t) + nC(t) + —j— — e ) h — 



dA w dt 2 ^ v ' \ 2 ^ v ' t J dh 
in which we now have set: 



W (n) =0, (3.14) 



d=2 + e. (3.15) 

The two RG functions can be obtained at one-loop order from a calculation of 
the 2-point function T^ 2 ); 

T (2) (p) = ^(p 2 + h) + [p 2 + \ (N - 1) h] ^ ^ + O(t) . (3.16) 

Applying the RG equation ( |3.1C| ) to r^ 2 ) and identifying the coefficients of p 2 
and n, we derive two equations which determine 3{t) and Q (t) at one-loop order 

3(t) =et- ~ 2 \ 2 + (t 3 , t 2 e) , (3.17a) 
CW = ^^t + 0(t 2 ,fe). (3.176) 



5.5 Discussion of the RG flow 

From the expression of 3(t) in equation (|3.17a| ) we immediately conclude: 
For d < 2 (e < 0), t = is an unstable IR fixed point, the IR instability being 
induced by the vanishing mass of the would-be Goldstone bosons. The spectrum 
of the theory thus is not given by perturbation theory and the perturbative as- 
sumption of spontaneous symmetry breaking at low temperature is inconsistent. 
As mentioned before, this result agrees with rigorous arguments. Note that since 
the model depends only on one coupling constant, t = is also a UV stable fixed 
point (the property of large momentum asymptotic freedom). Section 3J3 con- 



tains a short discussion of the physics in two dimensions for N > 2. The abelian 
case A^ = 2 is special and has to be discussed separately. 
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For d > 2, i.e. e > 0, t = is a stable IR fixed point, the O(N) symmetry 
is spontaneously broken at low temperature in zero field. The effective coupling 
constant, which determines the large distance behaviour, approaches the origin 
for all temperatures t < t c , t c being the first non-trivial zero of f3(t). Therefore 
the large distance properties of the model can be obtained from low temperature 
expansion and renormalization group, replacing the perturbative parameters by 
effective parameters obtained by solving the RG equations. 

The critical temperature. Finally we observe that, at least for e positive and 
small, and iV > 2, the RG function f3(t) has a non-trivial zero t c : 



2ne 



N 



+ O (e 2 ) (t c ) = , and (3' (t c ) = -e + O (e 2 ) . (3.18) 



Since t c is an unstable IR fixed point, it is by definition a critical temperature. 
Consequences of this property are studied below. Let us only immediately note 
that t c is also a UV fixed point, i.e. it governs the large momentum behaviour of 
the renormalized theory. The large momentum behaviour of correlation functions 
is not given by perturbation theory but by the fixed point. As a consequence 
the perturbative result that the theory cannot be rendered finite for d > 2 with 
a finite number of renormalization constants, cannot be trusted. 
We now discuss more precisely the solutions of the RG equations. 



3.4 Integration of the RG equations 

We first examine the implications of the RG equations for the large distance 
behaviour of correlation functions for d > 2 where t = is an IR fixed point. 
Equation (|3.10| ) can be solved as usual by the method of characteristics, i.e. by 
introducing a scaling parameter A and running parameters. It is here convenient 
to proceed somewhat differently by looking for a solution of the form 



r (n) ( Pi ; t, h, A) = C d (t)M - n (t)F^ ( Pi £(t), h/h (t)). 



(3.19) 



The ansatz ( |3.19|) solves the RG equations provided the unknown functions 
Mo(t), £(t) and ho(t) are chosen to be 



M (t) = exp 



1 



C(0 



2io P(t') 



f(t) = A-V/ S exp 



with then 



dt' 



1 



1 



\P{t') et> 



h (t) = tMo 1 (t)C d (t)A 



2-d 



dt' 



(3.20) 
(3.21) 

(3.22) 



Note that the function £(£) has in zero field the nature of a correlation length. 
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For the connected correlation functions the same analysis leads to: 

WW ( Pi ; t, H, A) = e d(n - 1} (t)M£{t)G< n) h/ho(t)) . (3.23) 

It is convenient to also introduce the function K(t) 



K(t) = M (t) [Amf' 2 /* = ! + 0(t). (3.24) 

Combining equation ( |3.19|) with dimensional analysis we can rewrite the scaling 
relation in an equivalent form 

T^( Pl ,t,h,A) ^ M^ n (t)[K(t)h] d / 2 

xrC») ( » tm^ (K_Y~ 2V2 1 11(3 25) 
1 ^[A-(t)fc]Va' M (t) UV 1 " 

Let us apply this result to the determination of the singularities near the coex- 
istence curve, i.e. at t fixed below the critical temperature when the magnetic 
field h goes to zero. 

The coexistence curve. The magnetization is given by 

M(t, h, A) = (a(x)) = A _e *-^-, (3.26) 

(r^ is the magnetic field dependent free energy). At one- loop order in a field 
one finds 

N — 1 . 1 f A d d q „ , 9N 

M = 1 A" e t / „ \ + O (t 2 ) . 

2 (27r) d 7 ? 2 + /i K J 

Thus from relation ( |3.25|) follows 
M(t, M = 1) = M (t) - [K(t)] d/2 /i (d " 2) / 2 r( ( 1 4 ~f/ 2 2) + O (h, h d -*) . 

This result shows that Mo(t) is the spontaneous magnetization and displays the 
singularity of the scaling equation of state (section |2l|) on the coexistence curve 
(H = 0, T < T c ) for N > 1, and in all dimensions d > 2. 

TTie equation of state in the critical domain. Let now instead use the scaling 
form (F7T9D 

M = A 2 "^^ = M (t)F^(h/h (t)). (3.27) 
Inversion of this relation yields the scaling form of the equation of state: 

h = Mi)/ {m) ■ (3 ' 28) 

28 



and the 1PI correlation functions can thus be written in terms of the magneti- 
zation as: 



r (B) (pu t, m, a) = r d (t) V(*)f (n) (nt (*) > J^j) • (3-29) 



The equations ( p.28| , |3.29[) are consistent with the equations ( |2.45| , |2.50| ): the ap 



pearance of two different functions £(t) and M (t) corresponds to the existence 
of two independent critical exponents u, f3 in the (4> 2 ) 2 field theory. They extend, 
in the large distance limit, the scaling form of correlation functions, valid in the 
critical region, to all temperatures below t c . There is however one important dif- 
ference between the RG equations of the (</> 2 ) theory and of the a- model: the 
((/> 2 ) theory depends on two coupling constants, the coefficient of <p 2 which plays 
the role of the temperature, and the coefficient of (</> 2 ) 2 which has no equivalent 

here. The correlation functions of the continuum (V> 2 ) 2 theory have the exact 
scaling form ( |3.29| ) only at the IR fixed point. In contrast, in the case of the 
ex-model, it has been possible to eliminate all corrections to scaling corresponding 
to irrelevant operators order by order in perturbation theory. We are therefore 
led to a remarkable conclusion: the correlation functions of the O(N) non-linear 
model are identical to the correlation functions of the (<^ 2 ) field theory at the IR 
fixed point. This conclusion is supported by the analysis of the scaling behaviour 
performed within the 1/N expansion (see equation ( f4.53|) ). 

Critical exponents. Let us now study more precisely what happens when t 
approaches t c (for N > 2). The function diverges as: 

~ A -1 (t c - f) 1 /* 9 '^) . (3.30) 
We conclude that the correlation length exponent v is given by 

v= ~m- (a31) 

For d close to 2 the exponent v thus behaves like: 

i/~l/e. (3.32) 

The function Mo(t) vanishes at t c : 

\nM (t) = ~j^\n(t c -t)+ const.. (3.33) 

This yields the exponent f3 and thus also r\ through the scaling relation (3 = 

V = C(t c )-e. (3.34) 
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A leading order we find: 

V = J^~ 2 + O (e 2 ) ■ (3.35) 

We finally note that the singularity of coming from the prefactor £ -d MQ~ n 
indeed agrees near £ c with the result of equation ( p.37|) . 

Consideration of operators with four derivatives allows also to calculate the 
exponent uj which characterizes leading corrections to scaling. One finds 

u = 4-d- 2e/(N - 2) + 0(e 2 ) . 

The nature of the correlation length £(£). The length scale £(£) is a cross-over 
scale between two different behaviours of correlation functions. For distances 
large compared to £(£), the behaviour of correlation functions is governed by 
the Goldstone modes (spin wave excitations) and can thus be deduced from the 
perturbative low temperature expansion. However when t approaches £ c , £(£) 
becomes large. There then exist distances large with respect to the microscopic 
scale but small with respect to £(£) in which correlation functions have a critical 
behaviour. In this situation we can construct continuum correlation functions 
consistent on all scales, the critical behaviour being also the large momentum 
behaviour of the renormalized field theory. 

General comment. From the consideration of the low temperature expansion 
we have been able to describe, for theories with a continuous symmetry, not only 
the complete structure of the low temperature phase, and this was expected, but 
also, in the non-abelian case, the critical behaviour near two dimensions . 

This result is somewhat surprising: Indeed perturbation theory is only sen- 
sitive to the local structure of the sphere S 2 = 1 while the restoration of sym- 
metry involves the sphere globally. This explains the peculiarity of the abelian 
case N = 2 because locally a circle cannot be distinguished from a non-compact 
straight line. For N > 2 the sphere has instead a local characteristic curvature. 
Still different regular compact manifolds may have the same local metric, and 
therefore the same perturbation theory. They all have the same low tempera- 
ture physics. However the previous results concerning the critical behaviour are 
physically relevant only if they are still valid when e is not infinitesimal and £ 
approaches £ c , a condition which cannot be checked directly. In particular the 
low temperature expansion misses in general terms decreasing like exp (const. /£) 
which may in some cases be essential for the physics. Therefore in section |4l] we 
will establish a direct relation between the O(N) a model and the (<^> 2 ) field 
theory to all orders in the large N expansion. This gives us some confidence 
that the previous considerations are valid for the iV-vector model at least for N 
sufficiently large. On the other hand the physics of N = 2 is not well reproduced. 
Cardy and Hamber have speculated about the RG flow for N close to 2 and di- 
mension d close to 2, incorporating phenomeno logically the Kosterlitz-Thouless 
transition in their analysis. 
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3.5 The dimension two 

Dimension two is of special interest from the particle physics point of view. 
The RG function j3{t) is then: 



(3(t) 




t 2 + (t 3 ) . 



(3.36) 



The non-linear cx-model for N > 2 is the simplest example of a so-called asymp- 
totically free field theory (UV free) since the first coefficient of the /3-function 
is negative, in contrast with the <p 4 field theory. Therefore the large momen- 
tum behaviour of correlation functions is entirely calculable from perturbation 
theory and RG arguments. There is, however a counterpart, the theory is IR 
unstable and thus, in zero field h, the spectrum of the theory is not perturba- 
tive. Contrary to perturbative indications, it consists of N massive degenerate 
states since the O(N) symmetry is not broken. Asymptotic freedom and the non- 
perturbative character of the spectrum are also properties of QCD, the theory 
of strong interactions, in four dimensions, . 
If we now define a function £(t) by: 



we can again integrate the RG equations and we find that £(t) is the correlation 
length in zero field. In addition we can use the explicit expression of the (3- 
function to calculate the correlation length or the physical mass for small t: 



However the exact value of the integration constant K, which gives the physical 
mass in the RG scale, can only be calculated by non-perturbative techniques. 

Finally the scaling forms ( |3.19| , |3~2~3| ) imply that the perturbative expansion at 
fixed magnetic field is valid, at low momenta or large distances, and for h/ho(t) 
large. 
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4 (4> 2 ) 2 Field Theory and Non-Linear a Model in the Large A Limit 

In the preceding sections we have derived universal properties of critical systems 
within the frameworks of the formal e = 4 — d and e = d — 2 expansions (at 
least for N > 2). It is therefore reassuring to verify, at least in some limiting 
case, that the results obtained in this way remain valid even when e is no longer 
infinitesimal. We show in this section that, in the case of the O(N) symmetric 
((/> 2 ) field theory, the same universal properties can also be derived at fixed 
dimension in the large iV limit, and more generally order by order in the large 
iV-expansion. We then examine the non-linear cr-model in the same limit. 

4-1 Introduction 
We again consider the partition function: 

Z = J [#(*)] exp [-5(0)], (4-1) 

where S(<p) is the O(N) symmetric action Q2.8Q (u = A 4 ~ d g): 

S{<t>) = j {^[a^(x)] 2 + ^ 2 (x) + ^[0 2 (x)] 2 |d d x. (4.2) 

A cut-off A, consistent with the symmetry, is implied. 

The solution of the model in the large A limit is based on a idea of mean field 
theory type: it can be expected that for A large the O(N) invariant quantities 
self-average and therefore have small fluctuations. Thus for example 

<0W(2/)> ~ (<P 2 (x))(<p 2 (y)). 

1 V — > oo 

This suggests to take (p 2 (x) as a dynamical variable. Technically, in the case 

of the (</> 2 ) theory, this can be achieved by using an identity similar to the 
Hubbard transformation: 



exp 



U> + 1 iff 



oc I dAexp(i-A 2 -^A-iA^), (4.3) 



where the integration contour is parallel to the imaginary axis. By introducing a 
field X(x) the identity can be used for each point x inside the functional integral 
( |4.1[ ). The new functional integral is then gaussian in <p and the integral over 
the field </> can be performed. The dependence on N of the partition function 
becomes explicit. Actually it is convenient to separate the components of <p into 
one component a, and A — 1 components tt, and integrate only over 7r (for T <T C 
it may even be convenient to integrate over only A — 2 components). For A large 
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the difference is negligible. To generate a correlation functions we also add a 
source H(x) to the action 



Z(H) = / [dX(x)} [d<r(x)]exp 



with: 



Sn (A, a) 



-S N (X,a) + J d d xH(x)a{x) 



+ 



(N-l) 



trln [-A + A(-)]- 



d d x 



(4.4) 



(4.5) 



X-field correlation functions. In this formalism it is natural to also calculate 
correlation functions involving the A-field. These have a simple interpretation in 
the initial (/>-field formalism. Indeed let us add a source j\ for A in the action 



(4.5). Then reintroducing the 0-field and integrating over A we recover instead 
of action (fO|), 

5(0) - (u/6)cf> 2 j x + (u/6)j 2 x - rj x . (4.6) 

Therefore j\ generates the 4> 2 correlation functions, up to a multiplicative factor 
and a translation of the connected 2-point function. 

4-2 Large N limit: the critical domain 

We now take the large iV limit at Nu fixed. With this condition Sn is of 
order iV and the functional integral can be calculated for N large by steepest 
descent. We expect a = 0(N 1 / 2 ), A = 0(1). We look for a uniform saddle point 
(a(x),X(x) space-independent), 

cr(x) = a , X(x) = X . 

with respect to a and A we obtain the saddle 



Differentiating then action (^ 
point equations: 



Act = 0. 



a 
~N 



Nu 



(X-r) + 



(27T) 



d d p 
p 2 + X 



0. 



(4.7a) 
(4.76) 



Remark. In the large iV limit the leading perturbative contributions come 
from chains of "bubble" diagrams of the form displayed in figure 1. These dia- 
grams form asymptotically a geometrical series which is summed by the algebraic 
techniques explained above. 







Fig. 1 Leading diagrams in the limit N — > oo. 
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The low temperature phase. Equation (|4.7a|) implies either a = or A = 0. 
In the low temperature phase cr, the average value of the field, does not vanish. 
Equation Q4.7fc| ) then yields: 



a 2 6 1 f A d d p 

: r -7KZ^ (4-* 



N Nu (2n) d J p 2 

Note that this equation has solutions only for d > 2. This is a manifestation 
of the Mermin-Wagner-Coleman theorem: in a system with only short range 
forces a continuous symmetry cannot be broken for d < 2, in the sense that the 
average a of the order parameter necessarily vanishes. Physically the would-be 
Goldstone modes are responsible for this property: being massless, as we know 
from general arguments and as the propagator in the r.h.s. of Q4.8|) confirms, 
they induce an IR instability for d < 2. 
Setting 

Nu 1 f A d d p 



we can rewrite equation 



a 2 



• , , ■ (4-9) 
r = r c +(u/6)T, (4.10) 



-r = (-r) 2f3 with (3= \ ■ (4.11) 



The expectation value of the field vanishes for r = r c , which therefore corresponds 
to the critical temperature. Moreover we find that for iV large the exponent f3 
remains classical, i.e. mean-field like, in all dimensions. 



The high temperature phase. Above T c , a vanishes. In expression ( [4.5|) we see 
that the a-propagator then becomes 

A a = ^—. (4.12) 
p z + A 

Therefore A 1 / 2 is at this order the physical mass, i.e. the inverse correlation 
length £ _1 of the field a 

m = C 1 = A 1/2 . (4.13) 

From equation ( f4.7fc| ) we can verify that dr/dX is positive. The minimum value 
of r, obtained for A = 0, is r c . Using equations ( fO| , |4.1Up in equation ( |4.7fcp we 
then find: 



6 + jv r a" v =JL _ 

u (27r) d J P 2 (p 2 + m2 ) m2 
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(i) For d > 4 the integral in ( }4.14j) has a limit for m = and therefore at 
leading order: 

m 2 = £ -2 ~ r and thus v = i , (4-15) 

which is the mean field result. 

(ii) For 2 < d < 4 instead, the integral behaves for m small like (setting 
d = 4-e): 



D 1 (m 2 ) = — ^ [ 2 . d 2 P = C(d)m- £ - a(d)A~ s + O (m 2 - £ A" 2 ) , 
(2n) d J P 2 (P 2 + m 2 ) K ' 

(4.16) 



(27T) 

with 



N * = t^mm <417o) 

where we have introduced for convenience the usual loop factor iV^. The con- 
stant a(d) which characterizes the leading correction in equation (|4.16 ), depends 
explicitly on the regularization, i.e. the way large momenta are cut. 

The leading contribution, for m —>■ 0, to the l.h.s. of equation fl4.14j ) now comes 
from the integral. Keeping only the leading term in (|4.16|) we obtain: 

m = ^- 1 -r 1 ^ 2 -^, (4.18) 

which shows that the exponent v is not classical: 

v=^— = -^—. (4.19) 
2-e d-2 K ' 

(iii) For d = 4 the l.h.s. is still dominated by the integral: 

^ / 2^ 1 f A d 4 p 1 , . A . . 

D\{m ) = - — tt / n , o ?rr ~ — rln(A/TO). 

The correlation length no longer has a power law behaviour but instead a mean- 
field behaviour modified by a logarithm. This is typical of a situation where the 
gaussian fixed point is stable, in the presence of a marginal operator. 

(iv) Examining equation ( |4.76p for a = and rf = 2we find that the correlation 
length becomes large only for r — ► — oo. This peculiar situation will be discussed 
in the framework of the non-linear cr-model. 



Finally, in the critical limit r = 0, A vanishes and thus from the form ( |4.12| ) of 
the cr-propagator we find that the critical exponent r\ remains classical for all d 

77 = d 4> = \(d-2). (4.20) 
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We verify that the exponents f3, u, r\ satisfy the scaling relation proven within the 
framework of the e-expansion 

(3 = vd<t> . 



Singular free energy and scaling equation of state. In a constant magnetic field 
H in the a direction, the free energy W(H)/Q per unit volume is given by 

3 V 1 M 

W(H)/Q = In Z/Q = — A 2 - —A - -Aa 2 + Ha - — tr(-A + A), 

2. 1 1 u 2 2 

where Q is the total space volume and A, a the saddle point values are given by 
equation (|4.7ft|) and the modified saddle point equation (|4.7aQ : 



Xa = H. (4.21) 
The thermo dynamical potential V(M) is the Legendre transform of W(H). First 

because partial derivatives of W with respect to A, a vanish as a consequence of 
the saddle point equations. It follows 

3 3r 1 TV 
V(M) = T(M) /D, = HM - W(H) /fi = A 2 + — A + - AM 2 + — tr(-A + A). 

As a property of the Legendre transform, the saddle point equation for A is now 
obtained by writing that the derivative of V vanishes. 

The term tr In can be evaluated for large A in terms of r c and the quantities 



defined in (|4.16|) . One finds 



trln[(A - A)A- X ] = J d d p ln[(p 2 + X)/p 2 } 



= _ 2 ^A d / 2 - ^A + ^A 2 A 4 " d + 0(A 1+d / 2 A- 2 ). 
d Nu 2 K 1 



The thermo dynamical potential becomes 



V(M) 4(i-I)A' + + -AM 2 - (4.22) 

v ' 2 \u* u) u 2 d ' v 1 

where we have defined 

u* = ^tt^tA 5 . (4.23) 
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Note that for A small the term proportional to A 2 is negligible with respect to 
the singular term \ d / 2 for d < 4. At leading order in the critical domain 

NC{d) 



V(M) 



-t\ + -XM 2 - 
2 2 



d 



-\ d " 



(4.24) 



where r has been defined in (|4.1C| ). 

The saddle point equation for A takes the simple form 

r + M 2 ~ NC(d)X d/2 - 1 = 0, 

and thus 

1 



A 



(r + M 2 ) 



2 /(d-2) 



_NC(d) 

It follows that the leading contribution, in the critical domain, to the thermody- 
namical potential is given by 

(d-2) 1 



V{M) 



2d 



(NC(d)) 



2 /(d-2) 



[ T + M 2 ) d ^ d -V 



(4.25) 



Various quantities can be derived from V(M), for example the equation of state 
by differentiating with respect to M. The resulting scaling equation of state is 

r)V 

' h M s f (t/M 2 ) , (4.26) 



H 



dM 



in which ho is a normalization constant, The exponent 6 is given by: 

. d + 2 
5 = ^2' 



(4.27) 
1, and the 



in agreement with the general scaling relation relation 5 = d/d^ 
function f(x) by: 

f(x) = (l + x) 2 / (d " 2) . (4.28) 

The asymptotic form of f(x) for x large implies 7 = 2/(d — 2) again in agreement 
with the scaling relation 7 = u(2 — 77). Taking into account the values of the 
critical exponents 7 and (3 it is then easy to verify that the function / satisfies 
all required properties like for example Griffith's analyticity (see section |2.6j ). In 
particular the equation of state can be cast into the parametric form: 

<T = R x ' % e, 
T = 3R(i-e 2 ), 

H = h R s/2 9 (3- 29 2 ) 2/{d - 2) . 




v - q 

Fig. 2 The "bubble" diagram Biy(p,m) 
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Leading corrections to scaling. The A 2 term yields the leading corrections to 
scaling. It is subleading by a power of r 



A 2 /A d/2 = 0(r (4 " d)/(d_2) ). 



We conclude 



to 



u= (4-d)/(d-2) 



uj = A - d. 



(4.29) 



We have identified the exponent u which governs the leading corrections to scal- 
ing. Note that for the special value u = u* this correction vanishes. 

Specific heat exponent. Amplitude ratios. Differentiating twice V(M) with 
respect to r we obtain the specific heat at fixed magnetization 



C 



H 



(d-2) (NC(d)f (d - 2) 
For M = we identify the specific exponent a 

4-d 



r + M 2 ) (4_d)/(d_2) . 



(4.30) 



a 



d-2' 



(4.31) 



which indeed is equal to 2 — du, as predicted by scaling laws. Among the ratio of 
amplitudes one can calculate for example and R c (for definitions see chapter 
28 of main reference) 



(R+) d 



4N T(3 - d/2) 
(d-2) 3 (47r) d / 2 



Rc 



4-d 
Jd^2~j'- 



(4.32) 



The A and (<p) 2 two-point functions. Differentiating twice the action (|4.5| ) with 
respect to X(x), then replacing the field X(x) by its expectation value m 2 , we 
find the A-propagator A\(p) above T c 



A,(P) = -jf 



6 

Nu 



+ B A (p : m) 



(4.33) 



where B\(p, m) is the bubble diagram of figure 2: 

"A Ad 



B A (p, m) 



1 



(27T) 



d d q 



(q 2 + m 2 ) (p — q) 2 + m 2 



(4.34) 



The A-propagator is negative because the A-field is imaginary. As noted in [|7T 
it is simply related to the <p 2 2-point function 



^ Y 1 1 + (Nu/6)B A (p,m) 



(4.35) 
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At zero momentum we recover the specific heat. The small m expansion of 
-E?a(0, m) can be derived from the expansion ( j4.16| ). One finds 



, . 1 f A d d q 

= ^(m 2 D 1 (m 2 )) = (d/2-l)C(d)m- £ - a(d)A" £ + • • • . (4.36) 

dm 1 m<A 

The singular part of the specific heat thus vanishes as m e , in agreement with 
equation flpq ) for M = 0. 

In the critical theory (m = at this order) for 2 < d < 4 the denominator is 
also dominated at low momentum by the integral 



B A (p, 0) = -Lf / 2 . d q = b(e)p- s - a(d)A~ s + O (A~ 2 p 2 ~ £ ) , 

(4.37) 



(2tt 

where 



and thus: 



sm(7ra/2) T(a — 1) 



Aa(p) ~ -l^P e - ( 4 -39) 
p^o Nb(e) 



We again verify consistency with scaling relations. In particular we note that in 
the large N limit the dimension [A] of the field A is 

[\] = l(d + e) = 2, (4.40) 

a result important for the 1/N perturbation theory. 
Remarks. 

(i) For d = 4 the behaviour of the propagator is still dominated by the integral 
which has a logarithmic behaviour oc l/ln(A/p). 

(ii) Note therefore that for d < 4 the contributions generated by the term 
proportional to A 2 (x) in (|4.5| ) always are negligible in the critical domain. 

4-3 RG functions and leading corrections to scaling 

The RG functions. For a more detailed verification of the consistency of the 
large iV results with the RG framework, we now calculate RG functions at leading 
order. One first easily verifies that, at leading order for A large, m solution of 
equation fl4.14j ) satisfies 

. dm „ T , , u 2 dm 
A_ +J V e a(d)A-» Y -=0, 
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where the constant a(e) has been defined in ( }4.16| ). It depends on the cut-off 
procedure but for e = 4 — d small satisfies 



a(e) ~ 1/(8ttV). (4.41) 

We then set (equation (|4.23| )): 



u = gA £ , g* = u*A~ £ = 6/(Na) . (4.42) 

In the new variables A, g, r we obtain an equation which expresses that m is RG 
invariant 

A dA +(3 ( g) dg~~ mi - 9)T d^) m(T ' 9 ' A) = ° ' (443) 
with 

P{g) = -eg(l-g/g*), (4.44) 
v~\g) =2 + V2 (g) = 2- eg/g*. (4.45) 

When a(d) is positive (but this not true for all regularizations, see the discussion 
below), one finds an IR fixed point g*, as well as exponents u = e, and v~ x = 
d — 2, in agreement with equations ( |4.29| , f4~7T9| ). In the framework of the e- 



expansion u is associated with the leading corrections to scaling. In the large iV 
limit to remains smaller than two for e < 2, and this extends the property to all 
dimensions 2 < d < 4. 

Finally, applying the RG equations to the propagator ( |4.12|) , we find 



V (g) = , (4.46) 

a result consistent with the value ( |4 . 2 0|) found for r\. 

Leading corrections to scaling. From the general RG analysis we expect the 
leading corrections to scaling to vanish for u = u* . This property has already 
been verified for the free energy. Let us now consider the correlation length or 
mass m given by equation (|4.14[) . If we keep the leading correction to the integral 
for m small (equation ( [4.16| )) we find 

1 - — + NC(d)m- £ + O (m 2 - £ A- 2 ) = -^r , (4.47) 

u u* m 2 

where equation ( |4.42|) has been used. We see that the leading correction again 
vanishes for u = u* . Actually all correction terms suppressed by powers of order 
e for d — > 4 vanish simultaneously as expected from the RG analysis of the </> 4 
field theory. Moreover one verifies that the leading correction is proportional to 
(u — u*)t £ /( 2 ~ £ \ which leads to uju = e/(2 — e), in agreement with equations 
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In the same way if we keep the leading correction to the A-propagator in the 
critical theory (equation (|4.37|) ) we find: 



Aa(p) = ~ 



6 6 j/ \ _ 



(4.48) 



where terms of order A 2 and 1/N have been neglected. The leading corrections 
to scaling again exactly cancel for u = u* as expected. 

Discussion. 

(i) One can show that a perturbation due to irrelevant operators is equivalent, 
at leading order in the critical region, to a modification of the (</> 2 ) 2 coupling. 
This can be explicitly verified here. The amplitude of the leading correction to 
scaling has been found to be proportional to 6/Nu — a(d)A~ £ where the value 
of a(d) depends on the cut-off procedure and thus of contributions of irrelevant 
operators. Let us call u' the (4> 2 ) 2 coupling constant in another scheme where a is 
replaced by a'. Identifying the leading correction to scaling we find the relation: 

6A e 6A e 
Nu- a{d) = N^- a ' {d) ' 



homographic relation which is consistent with the special form ( |4.44j ) of the 
/3-function. 

(ii) The sign ofa(d) . It is generally assumed that a(d) > 0. This is indeed what 
one finds in the simplest regularization schemes, like the simplest Pauli-Villars's 
regularization where a(d) is positive in all dimensions 2 < d < 4. Moreover a(d) 
is always positive near four dimensions where it diverges like 

a(d) ' 



8ir 2 e' 

Then there exists an IR fixed point, non-trivial zero of the /^-function. For this 
value u* the leading corrections to scaling vanish. 

However for d fixed, d < 4, this is not a universal feature. For example in the 
case of simple lattice regularizations it has been shown that in d = 3 the sign is 
arbitrary. 

However, if a(d) is negative, the RG method for large N (at least in the per- 
turbative framework) is confronted with a serious difficulty. Indeed the coupling 
flows in the IR limit to large values where the large N expansion is no longer 
reliable. It is not known whether this signals a real physical problem, or is just 
an artifact of the large N limit. 

Another way of stating the problem is to examine directly the relation between 
bare and renormalized coupling constant. Calling g r m 4 ~ d the renormalized 4- 
point function at zero momentum, we find 

™ 4 9r= l + A 4 -VVB A (0,m)/6' (449) 
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In the limit m <C A the relation can be written 



1_ = (d-2)NC(d) /mx4-d /i _ Na(d)_\ 

9r 12 \AJ \g 6 ; " 1 ' J 

We see that when a(d) < the renormalized IR fixed point value cannot be 
reached by varying g > for any finite value of m/A. In the same way leading 
corrections to scaling can no longer be cancelled. 

4-4 Small coupling constant and large momentum expansions for d < 4 

Section |4.6| is devoted to a systematic discussion of the 1/N expansion. How- 
ever the 1/N correction to the two-point function will help us to immediately 
understand the problem of the massless field theory for d < 4. 

We have seen that, in the framework at the 1/N expansion, we can calculate 
at fixed dimension d < 4 in the critical limit (T = T c , m 2 = 0). This implies 
that the terms of the 1/N expansion cannot be expanded in a power series of 
the coupling constant, at least with integer powers. Note that since the gaussian 
fixed point is an UV fixed point, the small coupling expansion is also a large 
momentum expansion. To understand the phenomenon we consider the (o~a) 
correlation function at order 1/N. At this order only one diagram contributes 
(figure 3), containing two X 2 a vertices. After mass renormalization and in the 
large cut-off limit we find: 

r( 2)/ ) = 2, 2 [ d d q f 1 _l\, n (J_ 

a*{P) V + N ( 2n )d J (<o/Nu) +b(e)q- s \(p + <l) 2 Q 2 J \N 2 

(4.51) 

An analytic study of the integral reveals that it has an expansion of the form 

a k u k p 2 - k£ + [3kU i2+2k)/£ p- 2k . (4.52) 

fe>i 

The coefficients oek, ftk can be obtained by performing a Mellin transformation 
over u on the integral. Indeed if a function f(u) behaves like u for u small, then 
the Mellin transform M(s) 



n OO 

M(s) = / duu^-'fiu), 
Jo 



has a pole at s = t. Applying the transformation to the integral, and inverting 
q and u integrations we have to calculate the integral 

u' 1 - 3 N ( Nb(e)q- £ \ 1 ~ s it 

du 



o (6/Nu) + b(e)q e 6 \ 6 / sin7rs 
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Fig. 3 The diagram contributing to T a a at order 1/N 



Then the value of the remaining q integral follows from the generic result ( |4.78| ) . 

The terms with integer powers of u correspond to the formal perturbative 
expansion where each integral is calculated for e small enough, has poles at 
e = (21 + 2)/k for which the corresponding power of p 2 is — /, i.e. an integer. 
One verifies that fti has a pole at the same value of e and that the singular 
contributions cancel in the sum. For these dimensions logarithms of u appear in 
the expansion. 



4-5 The non-linear a -model in the large N limit 

We have noticed that the term proportional to f d d x A 2 (a;), which has dimension 
4 — d for large N in all dimensions, is irrelevant in the critical domain for d < 4 
and can thus be omitted at leading order (this also applies to d = 4 where it 
is marginal but yields only logarithmic corrections). Actually the constant part 
in the inverse propagator as written in equation (|4.48|) plays the role of a large 
momentum cut-off. Let us thus consider the action QL67| ) without the A 2 term. 
If we then work backwards, reintroduce the initial field <p and integrate over X(x) 
we find 



[d<f>(x)] 5 



4> 2 {x) 



6 / 2 
— [m 

u v 



■) 



exp 



\{d^{x)f d d x 



(4.53) 



Under this form we recognize the partition function of the O(N) symmetric non- 
linear a-model in an unconventional normalization. We have therefore discovered 
a remarkable correspondence: to all orders in an 1/N expansion the renormalized 



non-linear a-model is identical to the renormalized 
fixed point. 



; ) field theory at the IR 



The large N limit. In order to more explicitly show the correspondence be- 
tween the set of parameters used in the two models, let us directly solve the 
c-model in the large N limit. We rewrite the partition function: 



[d(t>(x)d\(x)\ exp [S((f>, A)] 



(4.54) 



with: 



S(ct>, A) 



1 

2t 



d d x 



(d^f + A (</> 2 - 1) 



(4.55) 
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Integrating, as we did in section |4.1| , over N — 1 components of <f> and calling a 
the remaining component, we obtain: 



Z = J [da(x)d\(x)] exp [— Sn(ct, A)] 



(4.56) 



with: 



Sn (M) = i J [(^a) 2 + (a 2 (x)-l)A(x)] d d x + ~(iV - 1) tr In [-A + A(-)] • 

(4.57) 

The large iV limit is here taken at tN fixed. The saddle point equations, analo- 
gous to equations (|4.7|), are: 



m a = 



r2 = 1 . (N~l)t 



(2tt)< 



d d p 



pi _j_ m z 



(4.58a) 
(4.586) 



where we have set (X(x)) = m 2 . At low temperature a is different from zero and 
thus m, which is the mass of the 7r-field, vanishes. Equation (|4.586|) gives the 
spontaneous magnetization: 



1 - 



(N-l)t f A d d p 
(2n) d J ~f 



(4.59) 



Setting 



1 (N-l) f A d d p 



(2tt)« 



P 



we can write equation (|4.59|) 



(4.60) 



a 2 = l-t/t c . 



(4.61) 



Thus t c is the critical temperature where a vanishes. 

Above t c , a instead vanishes and m, which is now the common mass of the n- 
and cr-field, is for 2 < d < 4 given by: 



d d p 



1 1 _ md - 2 (^~ 1 ) 
t c t ' (2yr) d J p 2 (p 2 + 



- + (m 2 A d - 4 ) 
1 V 7 



(4.62) 



We recover the scaling form of the correlation length £ = 1/m. From the equa- 
tions ( [4.61| , |4~62D , we can also derive the RG functions at leading order for iV 
large: 



P(t) = et - ^-t 2 
Ztv 



Ztv 



(4.63) 
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It is also easy to calculate the thermo dynamical potential, Legendre transform 
of W(H) =t\nZ{H): 



V(M) = T(M)/0 = ^77^y(M 2 - 1 + t/t c ) d ^ d ~ 2 \ (4.64) 

2d (NC(d)) /( ' 

a result which extends equation ( f4.25| ) to all temperatures below t c . The cal- 
culation of other physical quantities and the expansion in 1/N follow from the 



considerations of previous sections and section [O. 



Two dimensions and the question of Borel summability. For d = 2 the critical 
temperature vanishes and the parameter m has the form: 



m ~ 



Ae" 2 ^^, (4.65) 



in agreement with the RG predictions. Note that the field 2-point function takes 
in the large iV-limit the form: 

T^(p)=p 2 + m 2 . (4.66) 

The mass term vanishes to all orders in the expansion in powers of the coupling 
constant t, preventing any perturbative calculation of the mass of the field. The 
perturbation series is trivially not Borel summable. Most likely this property is 
also true for the model at finite N. On the other hand if we break the O(N) 
symmetry by a magnetic field, adding a term ha to the action, the physical mass 
becomes calculable in perturbation theory. 



Corrections to scaling and the dimension four. In equation Q4.62|) we have 



neglected corrections to scaling. If we take into account the leading correction 
we get instead: 

m 2 (C(d)m d - 4 - a(d)A d - 4 ) xt-t e , 

where a(d), as we have already explained, is a constant which explicitly depends 
on the cut-off procedure and can thus be varied by changing contributions of 



irrelevant operators. By comparing with the results of section |4.3| , we discover 
that, although the non-linear a-model superficially depends on one parameter less 
than the corresponding (j) 4 field theory, actually this parameter is hidden in the 
cut-off function. This remark becomes important in the four dimensional limit 
where most leading contributions come from the leading corrections to scaling. 



For example for d = 4 equation (|4.62|) takes a different form, the dominant term 
in the r.h.s. is proportional to m 2 lnm. We recognize in the factor lnm the 
effective </> 4 coupling at mass scale m. Beyond the 1/N expansion, to describe 
with perturbation theory and renormalization group the physics of the non-linear 
a model it is necessary to introduce the operator J d d x\ 2 (x), which irrelevant 
for d < 4, becomes marginal, and to return to the </> 4 field theory. 
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4-6 The 1/N -expansion: an alternative field theory 

Preliminary remarks. Power counting. Higher order terms in the steepest 



descent calculation of the functional integral (|4.4|) generate a systematic 1/N 
expansion. Let us first slightly rewrite action (|L5p. We shift the field \(x) by 



its expectation value m (equation ( |4.13| )), \{x) 



m 



+ X(x) 



Sn (a, A) = i J d d x 
(N — 1 



(<V) 2 + mV + \{x)a 2 {x) - -\ 2 {x) - - (m 2 - r) \{x) 



u 



u 



+ 



trln [-A + m 2 + A( 



(4.67) 



We now analyze the terms in the action (|4.67| ) from the point of view of large iV 
power counting. The dimension of the field cr(x) is (d — 2)/2. From the critical 
behaviour ( |4.39| ) of the A-propagator we have deduced the canonical dimension 
[A] of the field A(x): 

2 [A] - e = d i.e. [A] = 2 . 

As noted above, A 2 has dimension 4 > d and is thus irrelevant. The interaction 
term J \(x)o~ 2 (x)d d x has dimension zero. It is easy to verify that the non-local 
interactions involving the A-field, coming from the expansion of the trln, have 
all also the canonical dimension zero: 



tr 



A(x) (-A + m 2 ) 



k [A] - 2k = . 



This power counting property has the following implication: In contrast with 
usual perturbation theory, the 1/N expansion generates only logarithmic correc- 
tions to the leading long distance behaviour for any fixed dimension d, 2 < d < 4. 
The situation is thus similar to the situation one encounters for the e-expansion 
(at the IR fixed point) and one expects to be able to calculate universal quanti- 
ties like critical exponents for example as power series in 1/N. However, because 
the interactions are non-local, the results of renormalization theory do not im- 
mediately apply. We now construct an alternative quasi-local field theory, for 
which the standard RG analysis is valid, and which reduces to the large iV field 
theory in some limit. 

An alternative field theory. To be able to use the standard results of renor- 
malization theory we reformulate the critical theory to deal with the non-local 
interactions. Neglecting corrections to scaling we start from the non-linear a- 
model in the form ( |4.55| ): 



[d\(x)} [d<t>(x)]exp l-S(<f),X)} 



1 



S{<f>, A) = f t I d d x 



+ A 



1) 



(4.68) 
(4.69) 
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The difficulty arises from the A-propagator, absent in the perturbative formu- 
lation, and generated by the large iV summation. We thus add to the action 



(|4.69| ) a term quadratic in A which at tree level of standard perturbation theory 
generates a A-propagator of the form ( |4.39|) . The modified action S v then is 



S v (<f>, A) 



d x , 

t 



(9, 



+ X 



1) 



rA(-a 2 



-e/2- 



(4.70) 



In the limit where the parameter v goes to infinity the coefficient of the additional 
term vanishes, and the initial action is recovered. 

We below consider only the critical theory. This means that the couplings of 
all relevant interactions will be set to their critical values. These interactions 
contain a term linear in A and a polynomial in <p 2 of degree depending on the 
dimension. Note that in some discrete dimensions some monomials become just 
renormalizable. We therefore work in generic dimensions. The quantities we 
shall calculate are regular in the dimension. The field theory with the action 
( 4.70 ) can be studied with standard field theory methods. The peculiar form 
of the A quadratic term, which is not strictly local, does not create a problem. 
Similar terms are encountered in statistical systems with long range forces. The 
simple consequence is that the A-field is not be renormalized because counter- 
terms are always local. 

It is convenient to rescale <fi i— ► <fiVt, A i— > v A: 



S v ((f>, X ) = \j ddx 



(9, 



A(— <9 2 ) e l 2 \ + relevant terms 



The renormalized critical action then reads 
1 



[S v 



d d x 



(d^f + v T Z v \<p 2 - A(-d 2 )" e/2 A + relevant terms 



(4-71) 

It follows that the RG equations for 1PI correlation functions of I A fields and n 
<p fields in the critical theory take the form: 



d d n 



p(Z,n) _ q _ 



(4.72) 



We can then calculate the RG functions as power series in 1/N. It is easy to verify 
that v 2 has to be taken of order 1/N. Therefore to generate a 1/N expansion 
one first has to sum the multiple insertions of the one-loop A two-point function, 
contributions which form a geometrical series. The A propagator then becomes 



4-d 



Aa(p) 



2p 
b(e)D(v) 



(4.73) 
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where we have defined 

D(v) = 2/b(e)+Nv 2 . 
The solution to the RG equations Ql72p can be written: 

r (z ' n) (rp,t;,A) = Z~ n//2 (r)r d ~ 2i_n(d ~ 2)/2 r (z ' n) (p, v(t), A), (4.74) 

with the usual definitions 

t^-=P(v(t)), t^^- = 7/(v(r)) . 
dr dr 

We are interested in the neighbourhood of the fixed point v 2 = oo. One verifies 
that the RG function 77(17) approaches the exponent rj obtained by direct calcula- 
tion, and the RG (3- function behaves like v 2 . The flow equation for the coupling 
constant becomes: 

dv 2 

t— = pv 2 , => v 2 {t) ~ r p . (4.75) 
dr 

We then note that to each power of the A field corresponds a power of v. It 
follows 

r( ,,n > (rj», v, A) oc ^( r ) r <*-2i-n(d-2+„) 

K T d-(2-p/2)l-n(d-2+ v ) ^ (4.76) 



To compare with the result ( 2.30 ) obtained from the perturbative renormalization 
group one has still to take into account that the functions r^' n ) defined here are 
obtained by an additional Legendre transformation with respect to the source of 
(p 2 . Therefore 

2- p/2 = d^2 = d-l/is. (4.77) 



Fig. 4 Diagram contributing to at order 1/N. 



I 




/ \ 
/ \ 
/ \ 
z 1 



Fig. 5 Diagram contributing to r^j A at order 1/N. 
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RG functions at order 1/N . Most calculations at order 1/N rely on the eval- 
uation of the generic integral 

1_ f d d q _ nd _ 2lM _ 2v Yjp + v- d/2)T(d/2 - g(g - u) _ 

(4.78) 



(2iv) d J (p + q) 2 vq 2 » P (47t) d / 2 T(iJ,)T(iy)T(d - \i - v 



For later purpose it is convenient to set: 

x = 2Nd = 4r(d-2) = 4sin(7r g /2)r(2- g ) 

1 b(e) T(d/2)T(2-d/2)T 2 (d/2-l) tiT(1 - e/2)T(2 - e/2) ' 1 ' J 

To compare with fixed dimension results note X\ ~ 2(4 — d) for d 4 and 
Xi ~ (d - 2) for d -»• 2. 

The calculation of the {4>4>) correlation function at order 1/N involves the 
evaluation of the diagram of figure 3. We want to determine the coefficient 
of p 2 lnA/p. Since we work at one-loop order we can instead replace the A 
propagator q~ e by q 2u and send the cut-off to infinity. We then use the result 
( 4.78Q with (j, = 1. In the limit 2v — > — e the integral has a pole. The residue 
of the pole yields the coefficient of p 2 In A and the finite part contains the p 2 lnp 
contribution 

Expressing that the function satisfies the RG equation we obtain the function 
n{v). 

The second RG function can be deduced from the divergent parts of the (<fi<fi\) 
function 



=v + Axv 3 D _1 (v) In A + A 2 v 5 D~ 2 (v) In A + finite 



with 



Ai = -r^rN d = -X, 
6(e) 

4N 

A 2 = -7o7^(d - 3)b(e)N d = —2N(d - 3)X X , 

where A\ and A 2 correspond to the diagrams of figures 4 and 5 respectively. 
Applying the RG equation one finds the relation at order 1/N 

P v *(v) = 2v 2 r](v) - 2A 1 v 4 D~ 1 (v) - 2A 2 v 6 D~ 2 (v). (4.80) 

We thus obtain 

9 

EV 

ri{v) = —X 1 D-\v) 1 (4.81) 

(3 v 2(v) = -^——X 1 D~ 1 (v) +4A^(1 - e)v 6 X 1 D~ 2 (v), (4.82) 
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where the first term in (3 V 2 comes from A\ and r\ and the second from A 2 . 
Extracting the large v 2 behaviour we find 



r ] = ^—^X 1+ 0(l/N 2 )., (4.83) 

4(3-e)(2- e ) 
P ~ N{A — e) Xl> °' 

\ = d - 2 + 2(3 ^ 2 ~ £) X 1 + 0(l/iV 2 ). (4.84) 



and thus 



4-7 Additional results 

The calculations beyond the order 1/N are rather technical. The reason is 
easy to understand: Because the effective field theory is renormalizable in all 
dimensions 2 < d < 4, the dimensional regularization, which is so useful in 
perturbative calculations, no longer works. Therefore either one keeps a true 
cut-off or one introduces more sophisticated regularization schemes. For details 
the reader is referred to the literature. 

Generic dimensions. The exponents 7 and r\ are known up to order 1/iV 2 and 
1/N 3 respectively in arbitrary dimensions but the expressions are too compli- 
cated to be reproduced here. The expansion of 7 up to order 1/N can be directly 
deduced from the results of the preceding sections: 

7 = 1=7/2 i 1 ~ W Xl ) +0 (jV 5 ) ' (486) 

The exponents u> and 6 = uju, governing the leading corrections to scaling, can 
also be calculated for example from the (A 2 AA) function: 

•-"-sM'-^MiW- (487) 

Note that the exponents are regular functions of e up to e = 2 and free of 
renormalon singularities at e = 0. 

The equation of state and the spin-spin correlation function in zero field are 
also known at order 1/N, but since the expressions are complicated we refer the 
reader to the literature for details. 

Three dimensional results. Let us give the expansion of rj in three dimensions 
at the order presently available: 

Vi , m . V3 . n ( 1 

7]= N + JP + W + \W 
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with 




tt 2 + f r (1/2) + 




7 



2 



24 



+ 



04 



( 




Ntt 2 



N 2 TT 4 



Note that the 1/N expansion seems to be rapidly divergent and certainly a 
direct summation of these terms does not provide very good estimates of critical 
exponents in 3 dimensions for useful values of N. 

4-8 Dimension four: triviality, renormalons, Higgs mass 

A number of issues concerning the physics of the (4> 2 ) 2 theory in four dimensions 
can be addressed within the framework of the large N expansion. For simplicity 
reasons we consider here only the critical (i.e massless) theory. 

Triviality and UV renormalons. It is easy to verify that the renormalized 
coupling constant g T , defined as the value of the vertex (aaaa) at momenta of 
order fx <C A, is given by: 



We see that when the ratio fi/A goes to zero, the renormalized coupling constant 
vanishes, for that all g. This is the so-called triviality property. In the standard 
treatment of quantum field field, one usually insists in taking the infinite cut-off 
A limit. Here one then finds only a free field theory. Another way of formulating 
the problem is the following: it is impossible to construct in four dimensions a 
4> field theory consistent (in the sense of satisfying all usual physical require- 
ments) on all scales for non zero coupling. Of course in the logic of effective field 
theories this is no longer an issue. The triviality property just implies that the 
renormalized or effective charge is logarithmically small as indicated by equa- 
tions (t4.88| , f3~89l) . Note that if g is generic (not too small) and A/li large, g r is 
essentially independent of the initial coupling constant. Only if the bare coupling 
is small is the renormalized coupling an adjustable, but bounded, quantity. 

Let us now imagine that we work formally and, ignoring the problem, we 
express the leading contribution to the four-point function in terms of the renor- 
malized constant: 



yr 1 + ±NgB A {n) ' 
where B\(p) corresponds to the bubble diagram (figure 2) 

B A (p) ~ — ^ln(A/p) + const.. 




(4.89) 



9 9r 



1 + HA/p) 1 + £^g T Hli/p) 
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We then find that the function has a poie for 

^/(N 9l ) 



This pole corresponds to the Landau ghost for this theory which has g = as an 
IR fixed point. If we calculate contributions of higher orders, for example to the 
two-point function, this pole makes the loop integrals diverge. In an expansion in 
powers of g T , each term is instead calculable but one finds, after renormalization, 
UV contributions of the type 



The perturbative manifestation of the Landau ghost is the appearance of contri- 
butions to the perturbation series which are not Borel summable. By contrast the 
contributions due to the finite momentum region, which can be evaluated by a 
semiclassical analysis, are Borel summable, but invisible for N large. This effect 
is called UV renormalon effect. Note finally that this UV problem is independent 
of the mass of the field (p, that we have taken zero for simplicity reasons. 

IR renormalons. We now illustrate the problem of IR renormalons with the 
same example of the massless (4> 2 ) 2 theory (but now zero mass is essential), in 
four dimensions, in the large N limit. We calculate the contribution of the small 
momentum region to the mass renormalization, at cut-off A fixed. In the large 
iV limit the mass renormalization is then proportional to (see equation ([4.51]) ) 



It is easy to expand this expression in powers of the coupling constant g. The 
term of order k in the limit k — > oo behaves as (— l) fc (iV/967r 2 ) fc /c!. This contri- 
bution has the alternating sign of the semiclassical contribution. Note that more 
generally for N finite on finds (— /32/2) fc /c!. IR singularities are responsible for 
additional, Borel summable, contributions to the large order behaviour. 

In a theory asymptotically free for large momentum, clearly the roles of IR 
and UV singularities are interchanged. 

The mass of the a field in the phase of broken symmetry. The </> 4 theory is 
a piece of the Standard Model, and the field a then represents the Higgs field. 
With some reasonable assumptions it is possible to establish for finite iV a semi- 
quantitative bound on the Higgs mass. Let us examine here what happens for 
iV large. 

In the phase of broken symmetry the action, after translation of average values, 
includes a term proportional to a A and thus the propagators of the fields a and 
A are elements of a 2 x 2 matrix M: 
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where a = (a(x)). In four dimensions Ba is given by equation (|4.89| ). It is 
convenient to introduce a mass scale M, RG invariant, such that 



487T 2 



+ 8tt 2 B a ( P ) ~ln(M/p), 



and thus 



Moce' 



The mass of the field a at this order is a solution to the equation det M = 0. 



The mass m a solution to the equation is complex, because the particle a can 
decay into massless Goldstone bosons. At a fixed, the mass decreases when the 
cut-off increases or when the coupling constant goes to zero. Expressing that the 
mass must be smaller than the cut-off, one obtains an upper-bound on m a (but 
which slightly depends on the chosen regularization) . 

4-9 Finite size effects 

Another question can be studied in the large N limit, finite size effects. It is 
difficult to discuss all possible finite size effects because the results depend both 
on the geometry of the system and on the boundary conditions. In particular one 
must discuss separately boundary conditions depending whether they break or 
not translation invariance. In the first case new effects appear which are surface 
effects, and that we do not examine here. Note that the periodic conditions are 
not the only ones which preserve translation invariance. For systems which have a 
symmetry one can glue the boundaries after having made a group transformation. 
Thus here one could also choose antiperiodic conditions or more generally fields 
differing by a transformation of the O(N) group. 

Moreover if we are interested only in qualitative aspects we can limit ourselves 
to a simple geometry, in each direction the system having the same finite size 
L, all other sizes being infinite (but we thus exclude some questions concerning 
crossover regimes). Even so the number of different possible situations remains 
large, and we limit ourselves here to two examples. 

We consider the example of periodic boundary conditions in two cases: finite 
volume (the geometry of the hypercube or rather hypertorus) in this section, and 
QFT at finite temperature in next section. 

From the point of view of renormalization group, finite size effects, which only 
affect the IR domain, do not change UV divergences. The RG equations remain 
the same, only the solutions are modified by the appearance of new dimensional 
quantities. Thus if finite sizes are characterized by only one length L, solutions 
will be functions of an additional argument L/£ where £ is the correlation length. 



One finds 



p 2 \n(M/p) 



(16tt 2 /AT)(t 2 m l\n(iM/m a ) = (16n 2 /N)a 2 . 
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A property characteristic of a system of finite size is the quantification of 
momenta in Fourier space. For periodic conditions, if we call L the size du 
system in each direction, we have 

p M = 2-kti^/L , n M G Z . 

In particular, in a massless theory the zero mode p = now corresponds to 
an isolated pole of the propagator. This automatically leads to IR divergences 
in all dimensions. Therefore in equations ( [4.71) the solution u / no longer 
exists. This is not surprising: there are no phase transitions in a finite volume. 
Neglecting corrections to scaling laws we can then write equation (|4.58£|) : 

1 = (N — l)tL~ d — - — — — , (4.90) 

where the sums are cut by a cut-off A. 

To discuss the equation it is convenient to introduce the function A(s) (related 
to Jacobi's elliptic functions) 

+oo 

A(s) = e ~ Sn2 ■ ( 4 - 91 ) 



n= — oo 



Using Poisson's transformation it is easy to show 

A(s) = (7t/s) 1/2 A (n 2 /s) . (4.92) 

Using this definition, and introducing the critical temperature t c , one can write 
equation (|4.9q ) (for 2 < d < 4) 

= (N — l)L~ d / ds (e- sm2 A d (A-K 2 s/L 2 ) - L d (4Tvs)- d/2 ) . (4.93) 

t t c Jq V / 

Setting s i— > L 2 s and introducing the function F: 

ds [e- sz2 A d (47t 2 s) - (47ts)- d/2 ) , (4.94) 

we can rewrite the relation 

- - - = (N- l)L 2 - d F(mL). (4.95) 

t t c 

For \t — t c \ <^ A d ~ 2 we find a scaling form which is in agreement with the RG 
result, which predicts = d-2 + 0{1/N)): 

Lm(t, L) = L/((t, L) = f((t- t c )L l /») . 
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Here the length £ has the meaning of a correlation length only for £ < L. Since 
77 = 0, the magnetic susceptibility x m zero held instead is always given by 
X = t/m 2 . 

One verifies that for t > t c fixed, L — > 00 and thus mL — > 00 one recovers 
the infinite volume limit. On the contrary in the low temperature phase for 
t < t c fixed, L — ► 00, mL goes to zero. Thus the contribution of the zero mode 
dominates in the r.h.s. of equation (|4.90| ). Using the relation (|4.92| ) one then 
finds 

F(z) = \ + K(d) + 0{z 2 ), 



K{d) 



ds 



A d (4n 2 s) - 1 - (4tts) 



-d/2 



and thus 



X (L, t) = -L = ^.(i _ t /t c )L d - tL 2 K(d) + O (L 4 - d /(t - Q) . (4.96) 
m z JM — 1 



We see that the susceptibility diverges with the volume, an indication of the 
existence of a broken symmetry phase. 

Note finally that it is instructive to make a similar analysis for different bound- 
ary conditions which have no zero mode. 

For d = 2 the regime where finite size effects are observables corresponds to 
t ln(LA) = O(l), i.e. to a regime of low temperature. The zero mode dominates 
for tln(LA) 1, and the susceptibility is then given by 

X (t,L)~±L 2 [1 + Q(t ln(LA))] . 



4-10 Field theory at finite temperature 

Quantum field theory at finite temperature can be considered as a system 
which has a finite size in one direction. Indeed the partition function is given 
by tre~ LH , where H is the hamiltonian and L _1 the temperature. For a scalar 
field theory with euclidean lagrangian density C(<fi) this leads to the functional 
integral 







j exp 


-J dr J d d ~ l x £((/>) 







where the field <f> satisfies periodic boundary conditions only in one direction 

4>(t = 0, x) = 4>{t = L, x). 

Let us again consider, as an example, the non-linear a model. We find a finite 
size system, but the interpretation of parameters is different. The variable t now 
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represents the coupling constant of the QFT. Since L is the inverse temperature, 
the limit L —>■ oo corresponds to the limit of vanishing temperature. 

The saddle point equation (|4.586|) , in the symmetric phase a = 0, becomes 



1 = (]V — l)t 



1 



(27T) 



d-1 



1 



m 2 + k 2 + (2-Kn/Lf 



(4.97) 



On immediately verifies that the IR problem induced by the zero mode has the 
following consequences: since one integrates only over d — 1 dimensions, a phase 
transition is only possible for d > 3. Qualitatively at large distance the condition 
of finite temperature leads to a property of dimensional reduction d ^ d — 1. The 
large iV expansion is thus particularly well suited to the study of this problem 
which exhibits a crossover between two different dimensions. 

Again using Schwinger's representation of the propagator, integrating over k 
and introducing the function ( f4.91| ) we can rewrite equation (^.97| ): 



1 



N-l 



2-d 



(47r)( d " 1 )/ 5 



G(mL) 



G{z) 



dss -(d-i)/2 [ e -* a M(47T 2 s) - (4ns)- 1 / 2 



(4.98) 
(4.99) 



Here £l 



m 



has really the meaning of a correlation length. 



This equation has a scaling form for d < 4. The behaviour of the system then 
depends on the ratio between L and the correlation length ^ of the system at 
zero temperature. For t > t c fixed and L large (with respect to 1/A) we recover 
the zero temperature limit. For t — t c small we find a crossover between a regime 
of small and high temperature. In the regime t < t c fixed and L large, we have 
to examine the behaviour of G(z) for z small. 

At d = 3: 

G(z) = — 2 In ,2 + const. . 



Hence 



1 



77r 



oc x{L, t) oc L 2 exp 



4nL 



1 



(4.100) 



One finds that £l remains finite below t c for all non vanishing temperatures, 
and has when the coupling constant t goes to zero or L — > oo the exponential 
behaviour characteristic of the dimension two. 

For d = 4 the situation is different because a transition is possible in dimension 
d — 1 = 3. This is consistent with the existence of the quantity G(0) > which 
appears in the relation between coupling constant and temperature at the critical 
point: 

1 1 _ (JV-l)G(O) 1 
t~T c ~ (4tt) 3 / 2 1? ' 



(4.101) 
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For a coupling constant t which corresponds to a phase of broken symmetry at 
zero temperature (£ < t c ), one now finds a transition temperature L _1 oc y/t c — t. 
Studying more generally the saddle point equations one can derive all other 
properties of this system. 



4-11 Other methods. General vector field theories 

The large iV limit can be obtained by several other algebraic methods. Without 
being exhaustive, let us list a few. Schwinger-Dyson equations for N large lead 
to a self-consistent equation for the two-point function. From the point of view of 
stochastic quantization or critical dynamics the Langevin equation also becomes 
linear and self-consistent for N large. One replaces <fi 2 (x,t) by (</> 2 (x,i)) ((■) 
means noise average) at leading order. Finally a version of the Hartree-Fock 
approximation also yields the large iV result. 

General vector field theories. We now briefly explain how the algebraic method 
presented in section [4.1| can be generalized to actions which have a more compli- 
cated dependence in one or several vector fields. Again in a general O(N) sym- 
metric field theory the composite fields with small fluctuations are the scalars 
constructed from all vectors. The strategy is then to introduce pairs of fields 
and Lagrange multipliers for all independent O(N) invariant scalar products 
constructed from the many-component fields. 

Let us first take the example of one field <p and assume that the interaction is 
an arbitrary function of the only invariant </> 2 (x) 

S(<f>) = J d d x{i[a^(x)] 2 + F(^ 2 )}. (4.102) 

We then introduce two fields p(x) and X(x) and use the identity: 



oc J [dp(x) dX(x)] exp |- J d d x [f A 2 -p) + V{pj\ J . 

(4.103) 

In the special case in which V(p) is a quadratic function, the integral over p can 
be performed. In all cases, however, the identity Q4.103|) transforms the action 
into a quadratic form in </> and therefore the integration over <p can be performed 
and the dependence in iV becomes explicit. This method will be applied in 
section [7] to the study of multi-critical points and double scaling limit. 

If the action is an O(N) invariant function of two fields <fii and 4>2 the potential 
depends on the three scalar products (j)\-<p%, 4>\ and 4>\. Then three pairs of fields 
are required. 



exp 



-J d d xV((f> 2 
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5 Gross— Neveu and Gross— Neveu— Yukawa Models 

To illustrate the techniques developed in sections [| f|, we now discuss models 
with fermions exhibiting the phenomenon of chiral phase transition. Again we 
consider two different field theory models with the same symmetries, the Gross- 
Neveu (GN) and the Gross-Neveu- Yukawa (GNY) models. The GN model is 
renormalizable in two dimensions, and describes in perturbation theory only one 
phase, the symmetric phase. The GNY model is renormalizable in four dimen- 
sions and instead allows a perturbative analysis of the chiral phase transition. 
We now show that the physics of these models can indeed be studied by the 
same techniques as ferromagnetic systems, that is RG equations near two and 
four dimensions, and large iV expansion. 

5.1 The Gross-Neveu model 

The GN model is described in terms of a U (N) symmetric action for a set of 
N massless Dirac fermions {-i/? 1 , 

$ ■ 0?P + ±G ■ ip) 2 . 

The GN model has in even dimensions a discrete chiral symmetry: 

tp^'js^, tp^-tjj'js, (5.1) 

which prevents the addition of a fermion mass term while in odd dimensions a 
mass term breaks space parity. Actually the two symmetry operations can be 
written in a form 

X = { Xl , X2, . . . , X d } -> X = {- Xl , X 2 , . . . , X d }, | ^ ^ _^ hi , 

valid in all dimensions. 

This model illustrates the physics of spontaneous fermion mass generation 
and, in even dimensions, chiral symmetry breaking. It is renormalizable and 
asymptotically free in two dimensions. However, as in the case of the non-linear a 
model, the perturbative GN model describes only one phase. The main difference 
is that the role of the spontaneously broken and the explicitly symmetric phase 
are interchanged. Indeed it is always the massless phase which is unstable in low 
dimensions. 

Since the symmetry breaking mechanism is non-perturbative it will eventually 
be instructive to compare the GN model with a different model with the same 
symmetries: the Gross-Neveu- Yukawa model. 

RG equations near and in two dimensions. The GN model is renormalizable 
in two dimensions, and in perturbation theory describes only the massless sym- 
metric phase. Perturbative calculations in two dimensions can be made with 



S(ip,1>) = - / d d x 
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an IR cut-off of the form of a mass term Mtjjijj, which breaks softly the chiral 
symmetry. It is possible to use dimensional regularization in practical calcula- 
tions. Note that in two dimensions the symmetry group is really 0(2N), as one 
verifies after some relabelling of the fields. Therefore the {ijji^) 2 interaction is 
multiplicatively renormalized. It is convenient to introduce here a dimensionless 
coupling constant 

u = GA 2 ~ d . (5.2) 

As a function of the cut-off A the bare correlation functions satisfy the RG 
equations: 



3 3 n 3 

A ^T + P( u )lT ~ nVA u ) ~ Vm(u)M- 



rW(p j ;u,M ) A) = 0. (5.3) 



3A ' ^ x '3u 2 ,VK ' llviy ' 3M 
A direct calculation of the /^-function in d = 2 + e dimension yields: 

/?(„) = en - (N< - 2)£ + (N' - 2)£ + (N ' ~ V + O («•) , (5.4) 

Note that for d = 2 N = 2N . 

The special case N' = 2, for which the /^-function vanishes identically in two 
dimensions, corresponds to the Thirring model (because for N' = 2 (i/j^^t/)) 2 = 
-2(^V>) 2 ). The latter model is to the equivalent the sine-Gordon or the 0(2) 
vector model. 

Finally the field and mass RG functions are 

N'-l 2 (N'-l)(N'-2) 3 (N'-l)(N'-2)(N'-5) 4 ( . 

^ {U) = 32^ U + 128^ ' (5 ' 5) 

. * N'-l N'-l 2 (2A'-3)(iV'-l) o 

As in the case of the non-linear a model, the solution of the RG equations 
(|5.3| ) involves a length scale £ of the type of a correlation length which is a RG 
invariant 

r f u du' 

£~V) = A(w) oc Aexp - J — — 



(5.6) 



Two dimensions. For d = 2 the model is asymptotically free. In the chiral the- 
ory (M — 0) the spectrum, then, is non-perturbative, and many arguments lead 
to the conclusion that the chiral symmetry is always broken and a fermion mass 
generated. From the statistical point of view this corresponds to the existence 
of a gap in the spectrum of fermion excitation (as in a superfluid or supercon- 
ductor). All masses are proportional to the mass parameter A (it) which is a RG 
invariant. Its dependence in the coupling constant is given by equation (|5.6| ): 

A(u) oc Au 1 /^'- 2 ^ e" 2 -/^'- 2 )" (1 + 0(u)) . (5.7) 
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We see that the continuum limit, which is reached when the masses are small 
compared to the cut-off, corresponds to u — > 0. 

S'-matrix considerations have then led to the conjecture that, for iV finite, the 
spectrum is: 

m n = A(u) 2 ^ N ~ ^ sin ^ J , ra = l,2...<iV,JV>2, 

To each mass value corresponds a representation of the 0(2N) group. The nature 
of the representation leads to the conclusion that n odd corresponds to fermions 
and n even to bosons. 

This result is consistent with the spectrum for iV large evaluated by semiclas- 
sical methods. In particular the ratio of the masses of the fundamental fermion 
and the lowest lying boson is: 

^= 2 Hw^) = 2+0{1/N2y (5 ' 8) 

The large N limit will be recovered in section |5.4| . 

Note that the two first values of iV are special, the model iV = 2 is conjectured 
to be equivalent to two decoupled sine-Gordon models. 

Dimension d = 2 + e. As in the case of the u-model, asymptotic freedom 
implies the existence of a non-trivial UV fixed point u c , in 2 + e dimension 



27T 

u r = — — -e 1 



N' — 2 \ N' 

u c is also the critical coupling constant for the transition between a phase in which 
the chiral symmetry is spontaneously broken and a massless small u phase. 
At the fixed point one finds the correlation length exponent v: 

z," 1 = -(3'{u c ) = e - -^—^ + O (e 3 ) . (5.9) 
The fermion field dimension [ip] is: 

N' — 1 

2M = d - 1 + ^K) = l + e+ 2(Ar/ _ 2)2 £ 2 + O {e 3 ) . (5.10) 

The dimension of the composite field a = ipip is given by 

[a] = d- 1 - rj M (u c ) = 1 



N' — 2 

As for the cr-model the existence of a non-trivial UV fixed point implies that large 
momentum behaviour is not given by perturbation theory above two dimensions, 
and this explains why the perturbative result that the model cannot be defined 
in higher dimensions cannot be trusted. However, to investigate whether the e 
expansion makes sense beyond an infinitesimal neighbourhood of dimension two 
other methods are required, like the 1/N expansion which will be considered in 
section FT4I. 
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5.2 The Gross-Neveu-Yukawa model 

The Gross-Neveu-Yukawa (GNY) model has the same chiral and U(N) sym- 
metries as the GN model. The action is (e = 4 — d): 

S ^,a) = J d d x -$-(0+ <7A e/2 a) V + I (<V) 2 + ±m V + ^AV 4 

(5.11) 

where a is an additional scalar field, A the momentum cut-off, and g, A dimen- 
sionless "bare" i.e. effective coupling constants at large momentum scale A. 

The action still has a reflection symmetry, a transforming into —a when the 
fermions transform by ( |5.1j ). In contrast with the GN model, however, the chi- 
ral transition can here be discussed by perturbative methods. An analogous 
situation has already been encountered when comparing the (4> 2 ) 2 field theory 
with the non-linear a model. Even more, the GN model is renormalizable in 
dimension two and the GNY model in dimension four. 



The phase transition. Examining the action ( |5.11| ) we see that in the tree 
approximation when m 2 is negative the chiral symmetry is spontaneously broken. 
The a expectation value gives a mass to the fermions, a mechanism reminiscent 
of the Standard Model of weak-electromagnetic interactions: 



m^=g(a), (5.12) 



while the a mass then is: 



m v = ^2 m l- ( 5 - 13 ) 



A 
37 2 

As a result of interactions the transition value m 2 of the parameter m 2 will be 
modified. In what follows we set 

m 2 = m 2 c +t, (5.14) 

where the new parameter t, in the language of phase transitions, plays the role 
of the deviation from the critical temperature. 

To study the model beyond the tree approximation we now discuss RG equa- 
tions near four dimensions. 

5.3 RG equations near four dimensions 



The model (|5 . 1 1|) is trivial above four dimensions, renormalizable in four di- 
mensions and can thus be studied near dimension 4 by RG techniques. Five 
renormalization constants are required, corresponding to the two field renormal- 
izations, the a mass, and the two coupling constants. The RG equations thus 
involve five RG functions. The 1PI correlation functions T^ l,n \ for I ip and n a 
fields, then satisfy 

( A + ^ W + ^ " * lv * " ^ " Vmt ^) Til,n) = ' (5-15) 
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Fig. 6 One-loop diagrams: fermions are represented by solid lines. 

The RG functions. The RG functions at one-loop order involve the calculation 
of the diagrams of figure 6. One finds: 



I3x = -e\ + 



1 



8tt s 



-A 2 + ANXg 2 - 24Ng l 



2 2iV + 3 4 



(5.16) 
(5.17) 



Note that in these expressions for convenience we have set in the algebra of 7 
matrices tr 1 = 4 as in four dimensions. To extrapolate the results to other 
dimensions one has to replace everywhere N by N'/4, where N' = iVtrl is the 
total number of fermion degrees of freedom. 

Dimension four. In four dimensions the origin A = g 2 = is IR stable. Indeed 
the second equation implies that g goes to zero, and the first then that A also 
goes to zero. As a consequence if the bare coupling constants are generic, i.e. if 
the effective couplings at cut-off scale are of order 1, the effective couplings at 
scale /i C A go to zero and in a way asymptotically independent from the bare 
couplings. One finds 



where we have defined 
A* 



8n< 



(2N + 3) ln(A///) 



A(/i) 



2 \* 



8tt 2 A 
ln(A/p) 



48N 



(2N + 3) [(2N - 3) + V4N 2 + 132N + 9j 



(5.18) 



This result allows to use renormalized perturbation theory to calculation physical 
observables. For example we can evaluate the ratio between the masses of the 
scalar and fermion fields. It is then optimal to take for [i a value of order (a). A 
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remarkable consequence follows: the ratio ( |5.13|) of scalar and fermion masses is 
fixed 

K = K = iw 

ml 3g 2 ( 2 jV - 3) + y/4N 2 + 132iV + 9 

while in the classical limit it seems arbitrary. 

Of course if the bare couplings are "unnaturally" small the same will apply to 
the renormalized couplings at scale /i and the ratio will be modified. 

Dimension d = 4 — e. One then finds a non-trivial IR fixed point (we recall 
N' = JVtrl): 

16n 2 e 

g2 * = wn' x * = Sn2eX *- ( 5 - 2 °) 

The matrix of derivatives of the /3-functions has two eigenvalues u, u', 



u 1 =e, u 2 = eVN' 2 + 132N' + 36/(N' + 6), (5.21) 

and thus the fixed point is IR stable. The first eigenvalue is always the smallest. 
The field renormalization RG functions are at the same order: 

N' 1 

= I^ 2 ' K> = Y^ g2 - (5 ' 22) 

At the fixed point one finds 

1* = WTS' ^ = (iV + 6) ' (5 - 23) 

and thus the dimensions and d a of the fields 

3 N' + 4 3 . . 

^=2"2(iV^) £ ' d ° = 1 -WT-6 £ - (5 ' 24) 

The RG function r\ m corresponding to the mass operator is at one-loop order: 

A 



Vm 16tt 2 V ° 



and thus the exponent v: 



1 e~ N'e 5N' + 6 + VN' 2 + 132JV' + 36 , 
- =2 + ,„ = 2- 3 A»- — =2- E . (5.25) 

Finally we can evaluate the ratio of masses ( |5.13|) at the fixed point: 
m 2 a A* 8JV' 



rn 



2 Q„2 



3g 2 (jy - 6) + VN' 2 + 132N' + 36 



In <i = 4 and d = 4— £ the existence of an IR fixed point has the same consequence: 
If we assume that the a expectation value is much smaller than the cut-off and 
that the coupling constants are generic at the cut-off scale, then the ratio of 
fermion and scalar masses is fixed. 
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5.4 GNY and GN models in the large N limit 

We now show that the GN model plays with respect to the GNY model 
(|5.11|) the role the non- linear a-model plays with respect to the <f> field the- 
ory. For this purpose we start from the action (|5.11|) of the GNY model and 
integrate over N — 1 fermion fields. We also rescale for convenience A^~ d ^ 2 ga 
into cr, and then get the large N action: 



S N ip,a)= / d d x < (0 + a) ip + A 



d-4 



1 fa \2 , m 2 , A 4 

|2^(V) + 7^ 

(iV-l)trln + 



^9 2g 2 A\g 

(5.26) 

To take the large N limit we assume a finite and g 2 , A = 0(1/N). 

Let us call I^(cr) the action per unit volume for constant field cr(x) and van- 
ishing fermion fields 

V(a) = A rf " 4 (^- 2 a 2 + ^a 4 ) - tftrln 0f + a) 

The expectation value of a for N large is given by a gap equation: 

V '(a)A 4 - d = ^a + —y - iV'A 4 ^^ t -^X, = . (5.28) 
V ' g 2 6g 4 (27r) d J q 2 + a 2 K ' 

It is also useful to calculate the second derivative to check stability of the extrema 

m 2 A 2 Ar/ .4_d f A d d g a 2 — q 2 



(2n) d (q 2 + a 



g 2 2^ 4 j 

The solution a = is stable provided 

™ 2 . 1 /" A d d 9 



2^2 



V"(0) > ^ ^ > iV'A 4 "^ / 



Instead the non-trivial solution to the gap equation exists only for 



g 2 (2yr) d 7 q 2 

but then it is stable. We conclude that the critical temperature or critical bare 
mass is given by: 

™g_ A r/ A 4-d 1 f Add( l 

g 2 (27t) d J q 



2 -^' A "7X> ( 5 - 29 ) 
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which shows that the fermions favour the chiral transition. In particular when d 
approaches 2 we observe that m 2 — ► +00 which implies that the chiral symmetry 
is always broken in 2 dimensions. Using equation ( |5.29| ) and setting 

t = A d - 4 (m 2 -m 2 c )/g 2 , (5.30) 

we can write the equation for the non-trivial solution 

t + hd-t—y + N'^ / A _4^-— = 0. 

6g 4 (2Tv) d J q 2 (q 2 +a 2 ) 

We now expand the integral for a small (equation ( [116] )) 

D 1 (a 2 ) = 7 2r2 q 2 =C(d)a- £ -a(d)A- s + 0[^—). (5.31) 

(2n) d J q 2 (q 2 +a 2 ) \ A 2 J 

Keeping only the leading terms for t —>■ we obtain for d < 4 the scaling 
behaviour 

a~ (-t/iVC) 1/(d - 2) . (5.32) 

Since, at leading order, the fermion mass = a, it immediately follows that 
the exponent v is also given by: 

v~P~ l/(d-2) => r] a = A-d. (5.33) 

At leading order, for iV — > 00, v has the same value as in the non-linear cr-model. 

At leading order in the scaling limit the thermo dynamical (or effective) poten- 
tial V(a) then becomes 

V(a) = \ta 2 + (N'/d)C(d)\a\ d . (5.34) 

We note that, although in terms of the a-field the model has a simple Ising-like 
symmetry, the scaling equation of state for large iV is quite different. 

We read from the large N action that at this order 77^ = 0. 

Finally from the large iV action we can calculate the a-propagator at leading 
order. Quite generally, using the saddle point equation, one finds for the inverse 
cr-propagator in the massive phase: 

N' f A d d q 

+ 2(27p ( ^ + 4a2) J (q 2 + a 2 )[(p + q) 2 + a 2 Y ^ 

We see that in the scaling limit p, a — > 0, the integral yields the leading contri- 
bution. Neglecting corrections to scaling we find that the propagator vanishes 
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for p 2 = —4a 2 which is just the 'ipt/j threshold. Thus, in this limit, m a = 2m^ in 
all dimensions, a result consistent with d = 2 exact value. 
At the transition the propagator reduces to 

(5.36) 



a N'b{e)p d - 2 
with (equation ( |4.38Q ) 

K^- .' (5.37) 
sm(7ra/2) F(d — 1) 

The result is consistent with the value of r) a found above. 

Let us finally note that the behaviour of the propagator at the critical point, 
A cr (p) oc p 2 ~ d , implies for the field a the canonical dimension [a] in the large N 
expansion, for 2 < d < 4: 

[a] = 1 . (5.38) 

Corrections to scaling and the IR fixed point. The IR fixed point is determined 
by demanding the cancellation of the leading corrections to scaling. Let us 
thus consider the effective potential V(o~). The leading correction to scaling is 
proportional to 



A N'a(d) 



4lg* 



(a (e) ~ l/87r 2 £). Demanding the cancellation of the coefficient of a 2 , we obtain 
a relation between A and g 2 

4 _ ^* _ 4A*£7T 2 ^ 2 \ 

9 * ~ 6N'a(d) ~ 3N' ^ 

a result consistent with the results of the e-expansion. 

In the same way it is possible to calculate the leading correction to the a- 
propagator (|5.35|) . Demanding the cancellation of the leading correction we 



obtain 

P ^ + ±-y-±N'(p 2 + 4o- 2 )a(d) = 0. 

The coefficient of a 2 cancels from the previous relation and the cancellation the 
coefficient of p 2 yields 

2 2 167T 2 £ ^ , 2N 

in agreement with the ^-expansion for large. 
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The relation to the GN model for dimensions 2 < d < 4. We have seen that 
the terms (d^a) 2 and cr 4 of the large iV action which have a canonical dimension 
4, are irrelevant in the IR critical region for d < 4. We recognize a situation 
already encountered in the (<f> 2 ) 2 field theory in the large iV limit. In the scaling 
region it is possible to omit them and one then finds the action: 



d d x 



-^■(0 + a)ij + A d -^o- 2 



(5.39) 



The integral over the a field can explicitly be performed and yields the action of 
the GN model: 



S N $,il>) = - / d d x 



2m z v ' 



The GN and GNY models are thus equivalent for the large distance physics. In 
the GN model, in the large iV limit, the a particle appears as a ipip boundstate 
at threshold. 

Conversely, it would seem that the GN model depends on a smaller number 
of parameters than its renormalizable extension. Again this problem is only 
interesting in four dimensions where corrections to scaling, i.e. to free field theory, 
are important. However, if we examine the divergences of the term tr In (0 + a) in 
the effective action (|5.26|) relevant for the large iV limit, we find a local polynomial 
in a of the form: 



d 4 



x 



Aa 2 (x) + B (c^cr) + Ca 4 (x) 



Therefore the value of the determinant can be modified by a local polynomial of 
this form by changing the way the cut-off is implemented: additional parameters, 
as in the case of the non-linear cr-model, are hidden in the cut-off procedure. Near 
two dimensions these operators can be identified with (i/ji/j) 2 , [d M ('(/>?/>)] 2 , (^i/j) 4 . It 
is clear that by changing the cut-off procedure we change the amplitude of higher 
dimension operators. These bare operators in the IR limit have a component on 
all lower dimensional renormalized operators. 

Note finally that we could have added to the GNY model an explicit breaking 
term linear in the a field, which becomes a fermion mass term in the GN model, 
and which would have played the role of the magnetic field of the ferromagnets. 

5.5 The large N expansion 

Using the large N dimension of fields and power counting arguments one can 
then prove that the 1/N expansion is renormalizable with arguments quite similar 



to those presented in section O . 



Alternative theory. To prove that the large iV expansion is renormalizable one 



proceeds as in the case of the scalar theory in section 4.6. One starts from a 
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critical action with an additional term quadratic in a which generates the large 
N a-propagator already in perturbation theory 



S(iP,iP,a) = Jd d x -*jj(0 + a)ij+ 7^(-d 2 ) d/2 -V 



(5.40) 



The initial theory is recovered in the limit v — > oo. One then rescales a in va. 
The model is renormalizable without a field renormalization because divergences 
generate only local counter-terms 



S v {^,^a) = J d d x 
RG equations follow 



-Z^(0 + v v Z v a)^ + ia(-a 2 ) d/2 -V 



(5.41) 



d d n 



p(/,n) _ q _ 



(5.42) 



Again the large N expansion is obtained by first summing the bubble contribu- 
tions to the a-propagator. We define 



Then the large N a propagator reads 



(aa) 



b(e)D{v)p d - 2 ' 
The solution to the RG equations can be written: 

T^ n \rp,v,A) = Z- n / 2 (T)r d - l - n{d - 2 ^ 2 T^ n \p,v(T), A), 
with the usual definitions 



(5.43) 



(5.44) 



dr 



dr 



We are interested in the neighbourhood of the fixed point v 2 = oo. Then the 
RG function rj(v) approaches the exponent rj. The flow equation for the coupling 
constant becomes: 

d^ 2 2 2/ \ o 

T— — = pV , =>• V (T) ~ T y . 

dr 

We again note that a correlation function with / a fields becomes proportional 
to Therefore 



r ( ' ,n) (rp, v, A) oc r d -( 1 -p/2)'-^(rf-2+^)/2_ 



(5.45) 
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We conclude 



d a = Ud - 2 + rj a ) = 1 - \p <S> i] a = 4 - d - p . (5.46) 



RG functions at order 1/N . A new generic integral is useful here 

1 f d d q^ + d) = , d - 2u - 2v r(n + u- d/2)T(d/2 l)r(d/2 - v) 

(2n) d J (p + q) 2 »q 2 » (47r) d / 2 r(//)r(i/)r(d - fi - u + 1) 

(5-47) 

We first calculate the 1/N contribution to the fermion two-point function at the 
critical point (from a diagram similar to diagram 3) 

r (2)/ ^ w , 2 ^ 2 f A d<W + flO 



' 6( £ )D(i;)(27r) d J Q d - 2 (p + Q) 2 ' 

We need the coefficient of ^lnA/p. Since we work only at one-loop order we 
again replace the a propagator l/q d ~ 2 by l/q 2v , and send the cut-off to infinity. 
The residue of the pole at 2u = d — 2 gives the coefficient of the term ^ln A and 
the finite part the tflnp contribution. We find 

,(2) , . ■ • Ziv 2 „ (d- 2 



r ^> } = ¥ + WW) Nd w J #ln(A/p) ' (5 ' 48) 

where N^ is the loop factor ( |4. 17q|) . Expressing that the Oipip) function satisfies 
RG equations we immediately obtain the RG function rj^(v) 

^ = W) (A ir Xl ' (5 ' 49) 

where X\ is given by equation ( [4.79Q . We then calculate the function (crtptp) at 
order 1/N 

r i ^(p) = v + A 1 D- 1 (v)v 3 lnA 1 

with 

6(e) 

where A\ corresponds to the diagram of figure 4. The diagram of figure 5 vanishes 
because the a 3-point function vanishes for symmetry reasons. 
The /3-function follows 

(3 v2 (v) = 4(d ~ 1) ^ X 1 g- 1 (^), (5.50) 
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and thus 

_ 8(d-l)N d _ 4(rf-l) 
9 db{e)N' dN' 1 ' 

The exponents rj^ and rj a at order 1/N, and thus the corresponding dimensions 
of fields d^, d a follow 

n _(d-2) Xl _ (d-2f T(d-l) 

d N' d r s (d/2)F{2-d/2)N r 1 ' } 

2d^=d-l- 2 -^§. (5.52) 

For d = 4 — e we find rj^ ~ s/N' , result consistent with ( |5.23| ) for A" large. For 
d = 2 + e instead one finds rj^ ~ e 2 /2N', consistent with ( |5.10|) . The dimension 
d a of the field a is 

d a = \{d - 2 + ria) = 1 - ^rT^i + 0(1/N' 2 ). (5.53) 

A similar evaluation of the (o" 2 crcr) function allows to determine the exponent 
v to order 1/N 

( ,54) 

1/ aAr 

Actually all exponents are known to order 1/N 2 except r)^ which is known to 
order 1/A^ 3 . 

6 Other models with chiral fermions 



Let us for completeness shortly examine two other models with chiral fermions, in 
which large A^ techniques can be applied, massless QED and the U(N) massless 
Thirring model. 

6.1 Massless electrodynamics 

Let us give another example with a structure different from a Yukawa-type 
theory. We now consider a model of A" charged massless fermion fields i/j, %[), 
coupled through an abelian gauge field (massless QED): 

S V, AJ = J d d x [^F 2 u (x) - $(x) • (0 + if) ^(x)] . (6.1) 

This model possesses, in addition to the U(l) gauge invariance, a chiral U(N) x 
U(N) symmetry because the fermions are massless. Again the interesting ques- 
tion is whether the model exhibits in some dimension 2 < d < 4 a spontaneous 
breaking of the chiral symmetry. 
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Dimension d = 4 — e. In terms of the standard coupling constant a: 



e 2 



(6.2) 

the RG j3 function reads (taking tr 1 = 4 in the space of 7 matrices) : 

2N 2 N o A(22A + 9) 4 
P(a) = -ea + —a 2 + + 7 r-^a 6 \ « 4 

3lV 27T 2 1447T 3 

1 iV[|i|iV 2 + (M6C(3)-^)iV + 23]a 5 + 0(a 6 ). (6.3) 
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The model is free at low momentum in four dimensions. Therefore no phase tran- 
sition is expected, at least for e 2 small enough. A hypothetical phase transition 
would rely on the existence on non-trivial fixed points outside of the perturbative 
regime. 

In the perturbative framework the model provides an example of the famous 
triviality problem. For a generic effective coupling constant at cut-off scale 
(i.e. bare coupling), the effective coupling constant at scale /i <C A is given 
by 

aKN ~ 4n ~ 2Nln(A/fj,y 

This result can be used to bound N. 

In 4 — e dimension, one instead finds a non-trivial IR fixed point corresponding 
to a coupling constant: 

,2 <-m_2 



(N' = iVtrl) and correlation functions have a scaling behaviour at long distance. 
As we have discussed in the case of the </> 4 field theory, the effective coupling 
constant at large distance becomes close to the IR fixed point, except when the 
initial coupling constant is very small. 

The RG function associated with the field renormalization is also known at 
order a 3 but this is a non-physical quantity since gauge dependent 

. a AN + 3 2 4(W 2 + 54N + 27 3 ^ , 4 , 
^ = ^-^6^ a + 576^ a 

where the gauge is specified by a term (<9 M A M ) 2 /2£. 

6.2 The large N limit 

To evaluate correlation functions for large, one first integrates over the 
fermion fields and one obtains the effective action: 

S(&lM„) = J d d x[^F^(x)-NtT\n(0 + i4)]. (6.4) 
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The large iV limit is taken with e 2 N fixed. Therefore, at leading order, only 
£2(^1)5 the quadratic term in in the expansion of the fermion determinant, 
contributes. A short calculation yields 

S 2 (A tl ) = -N' I & d kA„{k) [k 2 5^ - k^K] A v {-k)K{k), 

d-2 (6 ' 5) 
with K(k) = ^ _ - [b(e)k d - 4 - a(d)A d - 4 } + O (A" 2 ) , 

where a(d) is a regularization-dependent constant. 

For d < 4 the leading term in the IR region comes from the integral. The 
behaviour at small momentum of the vector field is modified, which confirms the 
existence of a non-trivial IR fixed point. The fixed point is found by demand- 
ing cancellation of the leading corrections to scaling coming from F 2 V and the 
divergent part of the loop integral, 



2 _ 2(d-l) A 4 " d 
(d-2)a(d) N' 



However there is still no indication of chiral symmetry breaking. Power counting 
within the 1/N expansion confirms that the IR singularities have been eliminated, 
because the large N vector propagator is less singular than in perturbation the- 
ory. Of course this result is valid only for N large. Since the long range forces 
generated by the gauge coupling have not been totally eliminated the problem 
remains open for d not close to four, or for e 2 not very small and N finite. Some 
numerical simulations indeed suggest a chiral phase transition for d = 4 and 
d = 3, N < N c ~ 3. 

The exponents corresponding to the IR fixed point have been calculated up to 
order 1/A 2 . At order 1/N (Xi is defined by equation ( |4.79| )) 

(d- 1)^(4 
n * d(d - 2) 4JV > 

, V) . 4 _ d _ (d - 3)(d - ( 6 J ( _ d - + O (!/*■) . 

Finally note that in the d = 2 limit, the integral generates a contribution 
Ne 2 /irk 2 times the propagator of the free gauge field 

N 1 

N'K(k) 

v J d^2 27tk 2 

As a direct analysis of the d = 2 case confirms, this corresponds to a massive 
bound state, of mass squared Ne 2 /n. However, for generic values of the cou- 
pling constant, this mass is of the order of the cut-off A. Only when e is small 
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with respect to the microscopic scale, as one assumes in conventional renormal- 
ized perturbation theory, does this mass correspond in the continuum limit to a 
propagating particle. 

Two dimensions. As stated above we now assume that the dimensional quan- 
tity e 2 is small in the microscopic scale. The model is then a simple extension of 
the Schwinger model and can be exactly solved in the same way. For N = 1 the 
model exhibits the simplest example of a chiral anomaly, illustrates the prop- 
erty of confinement and spontaneous chiral symmetry breaking. For N > 1 the 
situation is more subtle. The neutral ifjtfj two-point function decays algebraically 

($(x) ■ 1>(x)$(0) ■ V>(0)) oc x 2/N ~ 2 , 

indicating the presence of a massless mode and (4>ip) = 0. Instead if we calculate 
the two-point function of the composite operator On(x) 

N 

On(x) = Y\_ipi(x)ipi(x), 
i=i 

we find 

(O N (x)O N {0)) oc const. . 

We have thus identified an operator which has a non-zero expectation value. 
As a consequence of the fermion antisymmetry, if we perform a transformation 
under the group U(N) x U(N) corresponding to matrices U+, U-, the operator 
is multiplied by det U+ / det U- . The operator thus is invariant under the group 
SU(N) x SU(N) x £7(1). Its non-vanishing expectation value is the sign of the 
spontaneous breaking of the remaining U(l) chiral group. 

6.3 The U(N) Thirring model 
We now consider the model 

S$, 1>) = ~J $ (0 + m ) i> - \gJpJ,*] , (6.6) 

where 

^=^7m"V'- (6-7) 

The special case N = 1 corresponds to the simple Thirring model. In two 
dimensions it is then equivalent to a free massless boson field theory (with mass 
term for fermions one obtains the sine-Gordon model). Both to bozonize the 
model in d = 2 and to study that large N properties one introduces a abelian 
gauge field coupled to the current J M 

\gJ^J^ i— > Al/2g + J M . (6.8) 
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One then finds massive QED without the F£ v term 

S(A^ V) = - J d 2 x + i4 + ™o) ip - Al/2g] . (6.9) 

If we integrate over the fermions, the fermion determinant generates a kinetic 
term for the gauge field. For mo = we are thus in situation very similar to 
massless QED, except that the gauge field is massive. 
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7 The O(N) vector model in the large N limit: multi-critical points 
and double scaling limit 

We now discuss the large N limit of the general N- vector models with one scalar 
field. To illustrate the method we study multi-critical points. Of particular 
interest are the subtleties involved in the stability of the phase structure at 
critical dimensions. 

Another issue involves the so-called double scaling limit. Statistical mechan- 
ical properties of random surfaces as well as randomly branched polymers can 
be analyzed within the framework of large iV expansion. In the same manner 
in which matrix models in their double scaling limit provide representations of 
dynamically triangulated random surfaces summed on different topologies, O(N) 
symmetric vector models represent discretized branched polymers in this limit, 
where iV — > oo and the coupling constant g — > g c in a correlated manner. The 
surfaces in the case of matrix models, and the randomly branched polymers in 
the case of vector models are classified by the different topologies of their Feyn- 
man graphs and thus by powers of 1/N. Though matrix theories attract most 
attention, a detailed understanding of these theories exists only for dimensions 
d < 1. On the other hand, in many cases, the O(N) vector models can be suc- 
cessfully studied also in dimensions d > 1, and thus, provide us with intuition for 
the search for a possible description of quantum field theory in terms of extended 
objects in four dimensions, which is a long lasting problem in elementary particle 
theory. 

The double scaling limit in O(N) vector quantum field theories reveals an 
interesting phase structure beyond iV — > oo limit. In particular, though the 
N —>■ oo multicritical structure of these models is generaly well understood, 
there are certain cases where it is still unclear which of the features survives at 
finite N, and to what extent. One such problem is the multicritical behavior 
of O(N) models at critical dimensions. Here, one finds that in the iV — > oo 
limit, there exists a non-trivial UV fixed point, scale invariance is spontaneously 
broken, and the one parameter family of ground states contains a massive vector 
and a massless bound state, a Goldstone boson-dilaton. However, since it is 
unclear whether this structure is likely to survive for finite N one would like to 
know whether it is possible to construct a local field theory of a massless dilaton 
via the double scaling limit, where all orders in 1/N contribute. The double 
scaling limit is viewed as the limit at which the attraction between the O(N) 
vector quanta reaches a value at g — *■ g c , at which a massless bound state is 
formed in the iV — > oo limit, while the mass of the vector particle stays finite. In 
this limit, powers of 1/N are compensated by IR singularities and thus all orders 
in 1/N contribute. 

In section the double scaling limit for simple integrals and quantum me- 
chanics is explained, introducing a formalism which will be useful for field theory 
examples. 



80 



In section [7.2| the special case of field theory in dimension two is discussed. 

In higher dimensions a new phenomenon arises: the possibility of a sponta- 
neous breaking of the O(N) symmetry of the model, associated to the Goldstone 
phenomenon. 

Before discussing a possible double scaling limit, the critical and multicritical 
points of the O(N) vector model are examined in section |71^. In particular, 
a certain sign ambiguity that appears in the expansion of the gap equation is 
noted, and related to the existence of the IR fixed point in dimensions 2 < d < 4 



discussed in section OL In section 7.4 we discuss the subtleties and conditions 



for the existence of an O(N) singlet massless bound state along with a small 
mass O(N) vector particle excitation. It is pointed out that the correct massless 
effective field theory is obtained after the massive O(N) scalar is integrated out. 
Section |7l| is devoted to the double scaling limit with a particular emphasis 
on this limit in theories at their critical dimensions. In section |7T^ the main 
conclusions are summarized. 

7.1 Double scaling limit: simple integrals and quantum mechanics 

We first discuss d = and d = 1 dimensions, dimensions in which the matrix 
models has equally been solved. We however introduce a general method, not 
required here, but useful in the general field theory examples. 

The zero dimensional example. Let us first consider the zero dimensional 
example. The partition function Z is given by 



c z 



J d^exp [-NV (</> 2 )] 



The simplest method for discussing the large N limit is of course to integrate 
over angular variables. Instead we introduce two new variables A, p and use the 
identity 

exp [-NV(4> 2 )] oc J dpdAexp {-N [§A (<p 2 - p) + V{p)] } . (7.1) 

The integral over A is really a Fourier representation of a 5-function and thus 
the contour of integration runs parallel to the imaginary axis. The identity ( [7.1| ) 
transforms the action into a quadratic form in (p. Hence the integration over (p 
can be performed and the dependence in iV becomes explicit 

e z oc J dpdAexp{-iV [-\\p + V{p) + |lnA]}. 

The large iV limit is obtained by steepest descent. The saddle point is given by 

V'{p) = \\, P = i/\. 
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The leading contribution to Z is proportional to N and obtained by replacing 
A, p by the saddle point value. The leading correction is obtained by expanding 
A, p around the saddle point and performing the gaussian integration. It involves 
the determinant D of the matrix M of second derivatives 

M=("H" 2 v ~^j , D = detM = - 1 1 (V"(p)/X 2 + 1 1 ). 

In the generic situation the resulting contribution to Z is — \ In D. However if the 
determinant D vanishes the leading order integral is no longer gaussian, at least 
for the degree of freedom which corresponds to the eigenvector with vanishing 
eigenvalue. The condition of vanishing of the determinant also implies that two 
solutions of the saddle point equation coincide and thus corresponds to a surface 
in the space of the coefficients of the potential V where the partition function is 
singular. 

To examine the corrections to the leading large iV behaviour it remains how- 
ever possible to integrate over one of the variables by steepest descent. At leading 
order this corresponds to solving the saddle point equation for one of the vari- 
ables, the other being fixed. Here it is convenient to eliminate A by the equation 
A = 1/p. One finds 



e z oc 



J dp exp [-N(V(p) - \ lnp) + 0(1)} 



In the leading term we obviously recover the result of the angular integration 
with p = (p 2 . For N large the leading contribution arises from the leading term 
in the expansion of W(p) = V(p) — ^ lnp near the saddle point: 

W{p)-W{p a )~^W<ri{p B ){p-p a ) n . 

The integer n characterizes the nature of the critical point. Adding relevant 
perturbations 5 k V of parameters v k to the critical potential 

SkV = v k {p- p s ) k , l<k<n-2 

(the term k = n — 1 can always be eliminated by a shift of p) we find the partition 
function at leading order for N large in the scaling region: 

e z ^ Uk ^ oc J dz exp ^-z n - ^ u k z k 



where z oc A 1 / n (p — p a ) and 

u k oc N 1 ~ k / n Vk 
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is held fixed 

Quantum mechanics. The method we have used above immediately generalizes 
to quantum mechanics, although a simpler method involves solving the radial 
Schrodinger equation. We consider the euclidean action 



S(<j>)=N / dt ±(d t <f>(t)) 2 + V(<f 



(7.2) 



Note the unusual field normalization, the factor iV in front of the action simpli- 
fying all expressions in the large N limit. 

To explore the large iV limit one has to take the scalar function (p 2 , which 
self-averages, as a dynamical variable. At each time t we thus perform the 
transformation ( |7.1|) . One introduces two paths p(t), X(t) and writes 



exp 



N / dtVU 2 ) 



oc J [dp{t) dX(t)\ exp j-JV J dt [\X {cf 2 - p) + V{p)] j . (7.3) 
The integral over the path <p(t) is then gaussian and can be performed. One finds 



[dp(t)dX(t)} exp [S N (X,p)} 



(7.4) 



with 



S N =N J dt {-\Xp + V{p)\ +itrln(-d t 2 + A(-))- (7.5) 



Again, in the large N limit the path integral can be calculated by steepest 
descent. The saddle points are constant paths solution of 



1 



y'( P ) = ±A, P = — 



duo 



1 



2tt J u 2 + X 2sfX 



(7.6) 



where to is the Fourier energy variable conjugated to t. Again a critical point is 
defined by the property that at least two solutions to the saddle point equations 
coalesce. This happens when the determinant of the matrix of first derivatives 
of the equations vanishes: 



det I V Y 



(7.7) 



The leading correction to the saddle point contribution is given by a gaussian 
integration. The result involves the determinant of the operator second derivative 
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of Sn- By Fourier transforming time the operator becomes a tensor product of 
2x2 matrices with determinant D(u>) 

D(u) = det ( V "^ J , 1 with B(u) = -L f dU ' 



2 



^B(lu) J v ' 2tx J (a/ 2 + X)[(u - u') 2 + A] 



Thus, the criticality condition is equivalent to D(0) = 0. When the criticality 
condition is satisfied, the leading correction is no longer given by steepest descent. 
Again, since at most one mode can be critical, we can integrate over one of the 
path by steepest descent, which means solving the saddle point equation for one 
function, the other being fixed. While the p equation remains local, the A is now 
non-local, involving the diagonal matrix element of the inverse of the differential 
operator — d 2 + X(t). We shall see in next section how this problem can be 
overcome in general. A special feature of quantum mechanics, however, is that 
the determinant can be calculated, after a simple change of variables. We set 

X(t) = s(t) +s 2 (t), (7.8) 

in such a way that the second order differential operator factorizes 

-d 2 t + X(t) = -(d t + s(t)) (d t - s(t)). (7.9) 

The determinant of a first order differential operator can be calculated by expand- 
ing formally in s. Only the first term survives but the coefficient is ambiguous 

trln(l -d t _1 s(-)) = -0(0) J dts(t). 

A more refined analysis, which involves boundary conditions, is required to de- 
termine the ambiguous value 6(0) of step function. Here one finds 

lndet(-d? + A(-)) = tr ln(-d t 2 + A(-)) = / dts(t). (7.10) 

The jacobian of the transformation ( [7.8| ) contributes at higher order in 1/N and 
can be neglected. Therefore the large iV action becomes 

S N = N J dt [-\(s + s 2 )p + V(p) + |s] 

= N Jdt [-lps 2 + ls(p + l) + V(p)]. 

We can now replace s by the solution of a local saddle point equation (or perform 
the gaussian integration, but neglect the determinant which is of higher order) : 



8Sn 



^ -sp+ |(p + l) = 0. 
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and find 



'TV 



N / dt 



(7.11) 



We recognize the action for the large iV potential at zero angular momentum in 
the radial coordinate p(t) = 4> 2 (t). Critical points then are characterized by the 
behaviour of the potential W(p) 



W(p) = V{p) + 



1 



near the saddle point p s 



W(p)-W(p s )~w( n \p s ) 



(p-PsT 



nl 



At critical points the ground state energy, after subtraction of the classical term 
which is linear in N, has a non-analytic contribution. To eliminate N from the 
action we set 

t ^ t7V (n-2)/(n+2) ) ^ _ ^ ^ N -2/(n+2) ^ 

We conclude that the leading correction to the energy levels is proportional to 
jy-(n-2)/(n+2) ]\j Q ^ e a i so that the scaling of time implies that higher order time 
derivatives would be irrelevant, an observation which can be used more directly 
to expand the determinant in local terms, and will be important in next section. 
If we add relevant corrections to the potential 

hV = v k (p - p s ) k , l<fc<n-2, 

the coefficients Vk must scale like 



7.2 The 2D V(<p 2 ) field theory in the double scaling limit 

In the first part we study the O(N) symmetric V(4> 2 ) field theory, where cf> is iV- 
component field, in the large iV limit in dimension two because phase transitions 
occur in higher dimensions, a problem which has to be considered separately. 
The action is: 

S(4>)=N J d 2 x{i[^(x)] 2 + F(^ 2 )}, (7.12) 

where an implicit cut-off A is always assumed below. Whenever the explicit 
dependence in the cut-off will be relevant we shall assume a Pauli-Villars's type 
regularization, i.e. the replacement in action ( [7.1 2| ) of —(pd 2 <p by 

-<f)d 2 D(-d 2 /A 2 )(j), (7.13) 
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where D(z) is a positive non- vanishing polynomial with -D(O) = 1. 

As before one introduces two fields p(x) and A(x) and uses the identity (|7.3| ). 
The large N action is then: 

S N =N J d 2 x [V(p) - \Xp] + §JVtrln(-A + A). (7.14) 

Again for N large we evaluate the integral by steepest descent. Since the saddle 
point value A is the </>-field mass squared, we set in general A = m 2 . With this 
notation the two equations for the saddle point m 2 , p s = (<fi 2 ) are: 

V'{ Ps ) = \m 2 1 (7.15a) 



where we have used a short-cut notation 

i r A d 2 k i r d 2 k 



(2tt) 2 J k 2 + m 2 (2tt) 2 J D{k 2 / K 2 )k 2 + m? 
For m <C A one finds 



^B^m 2 ). (7.16) 



BUm 2 ) = — ln(A/m) + — ln(8nK) + 0(m 2 /A 2 ), 
27T Arc 

where K is a regular ization dependent constant. 

As we have discussed in the case of quantum mechanics a critical point is 
characterized by the vanishing at zero momentum of the determinant of second 
derivatives of the action at the saddle point. The mass-matrix has then a zero 
eigenvalue which, in field theory, corresponds to the appearance of a new massless 
excitation other than <f>. In order to obtain the effective action for this scalar 
massless mode we must integrate over one of the fields. In the field theory 
case the resulting effective action can no longer be written in local form. To 
discuss the order of the critical point, however, we only need the action for space 
independent fields, and thus for example we can eliminate A using the A saddle 
point equation. 

The effective p potential W(p) then reads 

/Hp) f) 
dX'X f —B 1 (X'), (7.17) 

where at leading order for A large 

X(p) = 8nKA 2 e- 47rp . 
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The expression for the effective action in equation ( [7.1 7|) is correct for any d and 
will be used also in section [7.51. Here we have: 



W(p) = V(p) + KA 2 e~ 47Tp = V(p) + ±m 2 e -*<p-Ps) . 

A multicritical point is defined by the condition 

W(p)-W( Ps ) = 0((p- Ps r) (7.18). 

This yields the conditions: 

V {k \p s ) = i(-4 7 r) fc - 1 m 2 forl<fc<n-l. 

Note that the coefficients V^ k \p s ) are the coupling constants renormalized at 
leading order for N large. If V(p) is a polynomial of degree n — 1 (the minimal 
polynomial model) the multicritical condition in equation ( 7.18|) determines the 
critical values of renormalized coupling constants as well as p s 

When the fields are space-dependent it is simpler to eliminate p instead, be- 
cause the corresponding field equation: 

V(p{x)) = \X{x). (7.19) 

is local. This equation can be solved by expanding p{x) — p s in a power series in 
X(x) — m 2 : 

p[x) ~ ps = wkp7) (A(x) " m2) + ((A " m2)2) ' (7 ' 20) 

The resulting action for the field X(x) remains non-local but because, as we shall 
see, adding powers of A as well as adding derivatives make terms less relevant, 
only the few first terms of a local expansion of the effective action will be impor- 
tant. 

If in the local expansion of the determinant we keep only the two first terms 
we obtain an action containing at leading order a kinetic term proportional to 
(<9 M A) 2 and the interaction (A(x) — m 2 ) n : 



1 



S N (X) ~ N J d 2 x 



:.r) - m 2 ) n 



where the neglected terms are of order (A — m 2 ) n+1 , A<9 4 A, and X 2 d 2 X and 

S n = W^{ Ps )[2V"{p 8 )]- n = W^(p s )(-4nm 2 )- n . 

Moreover we note that together with the cut-off A, m now also acts as a cut-off 
in the local expansion. 
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To eliminate the A dependence in the action we have, as in the example of 
quantum mechanics, to rescale both the field A — m? and space: 

A(x) - m 2 = v / 487rm 2 iV" 1/2 ^(x) , x i-> A (n " 2)/4 x . (7.21) 

We find 

S N ((p)~ J d 2 X [lidtf) 2 + fan*"] . 

In the minimal model, where the polynomial V(p) has exactly degree n — 1, we 
find g n = 6(487r)( n - 2 )/ 2 m 2 . 

As anticipated we observe that derivatives and powers of f are affected by 
negative powers of N, justifying a local expansion. However we also note that 
the cut-offs (A or the mass m) are now also multiplied by ]V~( n-2 )/ 4 . Therefore 
the large A limit also becomes a large cut-off limit. 

Double scaling limit. The existence of a double scaling limit relies on the 
existence of IR singularities due to the massless or small mass bound state which 
can compensate the 1/N factors appearing in the large A perturbation theory. 

We now add to the action relevant perturbations: 

S k V = v k (p(x) - Ps ) k , l<k<n-2. 
proportional to J d 2 x(A — m 2 ) : 

6 k S N (X) = NSk J d 2 x(A-m 2 ) fc , 
where the coefficients Sk are functions of the coefficients v k . After the rescaling 



5 k S N ((p) = 




However, unlike quantum mechanics, it is not sufficient to scale the coefficients 
g k with the power j\K fc_n )/ 2 to obtain a finite scaling limit. Indeed perturbation 
theory is affected by UV divergences, and we have just noticed that the cut-off 
diverges with N. In two dimensions the nature of divergences is very simple: 
it is entirely due to the self-contractions of the interactions terms and only one 
divergent integral appears: 

where p is the small mass of the bound state, required as an IR cut-off to define 
perturbatively the double scaling limit. We can then extract the N dependence 

(</? 2 (x)) = — (n - 2) In A + 0(1). 
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Therefore the coefficients Sk have also to cancel these UV divergences, and thus 
have a logarithmic dependence in N superposed to the natural power obtained 
from power counting arguments. In general for any potential U(tp) 



U(<p) =: Ultp) : + 



E^i^ 2 > fc ( 

k=i v 



k ( d 



2k' 



:U(<p) 



where : U(<p) : is the potential from which self-contractions have been subtracted 
(it has been normal-ordered). For example for n = 3 

<p 3 (x) =: ip 3 (x) : +3(v9 2 )v?(x), 

and thus the double scaling limit is obtained with the behaviour 

1 



N 9i + T7T- lniV#3 fixed . 

167T 



For the example n = 4 



tfiiV 3 / 2 and Ng 2 + — IniV fixed . 

57T 



7.3 The V{4> 2 ) in the large N limit: phase transitions 

In higher dimensions something new happens: the possibility of phase transi- 
tions associated with spontaneous breaking the O(N) symmetry. In the first part 
we thus study the O(N) symmetric V(<p 2 ) field theory, in the large N limit to ex- 
plore the possible phase transitions and identify the corresponding multicritical 
points. The action is: 

S(</>) = N J d d x { \ [d^(x)f + V (<P 2 ) } , (7.22) 

where, as above (equations (|7.12j , [7.13|) ), an implicit cut-off A is always assumed 
below. 

The identity ( [7.3| ) transforms the action into a quadratic form in <p and there- 
fore the integration over <fi can be performed. It is convenient however here to 
integrate only over N — 1 components, to keep a component of the vector field, 
which we denote a, in the action. The large N action is then: 

S N = N J d d x [| (<V) 2 + V{p) + f A (a 2 - p)] + ±(iV - 1) trln(-A + A). 

(7.23) 

The saddle point equations: the O(N) critical point. Let us then write the 
saddle point equations for a general potential V. At high temperature a = 
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and A is the </>-field mass squared. We thus set in general A = m 2 . With this 
notation the three saddle point equations are: 



m 2 £7 = , (7.24a) 

V'(p) = l m 2 , (7.246) 

1 f A c\ d k 

a 2 = p - -^—j / 9 . (7.24c) 

P (2n) d J k 2 + m 2 K ' 

In the ordered phase a ^ and thus m vanishes. Equation ( [7.24cj) has a solution 
only for p > p c , 

"A A d 



(27r) d 



J -p-> a = Vp - Pc 



Equation ( |7.24£| ) which reduces to V'(p) = then yields the critical temperature. 
Setting V(p) = U(p) + \rp, we find 

r c = -2U'(p c ). 

To find the magnetization critical exponent (5 we need the relation between the 
r and p near the critical point. 

In the disordered phase, a = 0, equation (|7.24c|) relates p to the 0-field mass 
m. For ro< A, p approaches p c , and the relation becomes (equation (|4.16|) ): 

p-p c = -C(d)m d - 2 + a(d)m 2 A d - 4 + O (m d A" 2 ) + O (m 4 A d " 6 ) . (7.25) 

For 2 < d < 4 (the situation we shall assume below except when stated otherwise) 
the O (m d A -2 ) from the non-analytic part dominates the corrections to the 
leading part of this expression. For d = 4 instead 

p — p c = g^ni 2 (lnm/A + const.) , 

and for d > 4 the analytic contribution dominates and 

p-p c ~ a(d)m 2 A d " 4 . 

The constant C(d) is universal (equation (|4.176|) ). The constant a(d), which 
also appears in equation ( fll6| ) , instead depends on the cut-off procedure, and is 
given by 

°w = wrl w f 1 - nw)) ■ (7 - 26) 

Critical point. In a generic situation V"(p c ) = U"(p c ) does not vanish. We 
thus find in the low temperature phase 

t = r _ rc „ -2U"(p c )(p -p c ) => (7.27) 
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This is the case of an ordinary critical point. Stability implies V"(p c ) > so 
that t < 0. 

At high temperature, in the disordered phase, the </>-field mass m is given by 
2U'(p) + r = m 2 and thus, using ( |7.25|) , at leading order 



t ~ 2U"( Pc )C(d)m d - 2 , 

in agreement with the result of the normal critical point. Of course the simplest 
realization of this situation is to take V(p) quadratic, and we recover the (4> 2 ) 2 
field theory. 

The sign of the constant a(d). A comment concerning the non-universal con- 
stant a(d) defined in (|7.25|) is here in order because, while its absolute value is 
irrelevant, its sign plays a role in the discussion of multicritical points. Actu- 
ally the relevance of this sign to the RG properties of the large N limit of the 
simple (4> 2 ) 2 field theories has already mentioned (section |4.3| ). For the simplest 
Pauli-Villars's type regularization we have D(z) > 1 and thus a(d) is finite and 
positive in dimensions 2 < d < 4, but this clearly is not a universal feature. 

A new situation arises if we can adjust a parameter of the potential in such 
a way that U"(p c ) = 0. This can be achieved only if the potential V is at 
least cubic. We then expect a tricritical behavior. Higher critical points can 
be obtained when more derivatives vanish. We shall examine the general case 
though, from the point of view of real phase transitions, higher order critical 
points are not interesting because d > 2 for continuous symmetries and mean- 
field behavior is then obtained for d > 3. The analysis will however be useful in 
the study of double scaling limit. 

Assuming that the first non- vanishing derivative is U^ n \p c ) : we expand further 
equation (|7.24£|) . In the ordered low temperature phase we now find 

t = -J^yU^( Pc )(p-p c r-\ =* P=^rTy (7-28) 

which leads to the exponent f3 expected in the mean field approximation for such 
a multicritical point. We have in addition the condition £A n ) (p c ) > 0. 
In the high temperature phase instead 



m 2 = t + (-l) n_1 7 r 7 C/' (n) (p c )C n - 1 (d)m (n - 1)(d - 2) . (7.29) 

(n — 1)! 



For d > 2n/(n — 1) we find a simple mean field behavior, as expected since we 
are above the upper-critical dimension . 

For d < 2n/(n — 1) we find a peculiar phenomenon, the term in the r.h.s. is 
always dominant, but depending on the parity of n the equation has solutions 
for t > or t < 0. For n even, t is positive and we find 

m0C ^ V= (n-l)(d-2Y (7 - 30) 
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which is a non mean-field behavior below the critical dimension. However for 
n odd (this includes the tricritical point) t must be negative, in such a way 
that we have now two competing solutions at low temperature. We have to 
find out which one is stable. We shall verify below that only the ordered phase 
is stable, so that the correlation length of the 0-field in the high temperature 
phase remains always finite. Although these dimensions do not correspond to 
physical situations because d < 3 the result is peculiar and inconsistent with the 
e-expansion. 

For d = 2n/(n — 1) we find a mean field behavior without logarithmic correc- 
tions, provided one condition is met: 

J^Ty U^HpJC"- 1 (2n/(n - 1)) < 1 , C(3) = l/(4rr). (7.31) 

We examine, as an example, in more details the tricritical point below. We will 
see that the special point 

^—-U{n) {pc)C n-i {2n/{n _!)) = !, (7.32) 
(n — 1)! 

has several peculiarities. In what follows we call O c this special value of (p c ). 

Discussion. In the mean field approximation the function U(p) oc p n is not 
bounded from below for n odd, however p = is the minimum because by 
definition p > 0. Here instead we are in the situation where U(p) ~ (p — p c ) n 
but p c is positive. Thus this extremum of the potential is likely to be unstable 
for n odd. To check the global stability requires further work. The question is 
whether such multicritical points can be studied by the large iV limit method. 

Another point to notice concerns renormalization group: The n = 2 example 
is peculiar in the sense that the large iV limit exhibits a non-trivial IR fixed 
point. For higher values of n no coupling renormalization arises in the large N 
limit and the IR fixed point remains pseudo-gaussian. We are in a situation 
quite similar to usual perturbation theory, the f3 function can only be calculated 
perturbatively in 1/N and the IR fixed point is outside the perturbative regime. 

Local stability and the mass matrix. The matrix of the general second partial 
derivatives of the effective action is: 

(p 2 + m 2 a \ 

V"(p) -\ , (7.33) 

a -\ -±B A (p,m) J 

where B\(p, m) is defined in (|4.34| ). 

We are in position to study the local stability of the critical points. Since 
the integration contour for A = m 2 should be parallel to the imaginary axis, a 
necessary condition for stability is that the determinant remains negative. 
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The disordered phase. Then a = and thus we have only to study the 2x2 
matrix M of the p, m 2 subspace. Its determinant must remain negative which 
implies 

det M < 2V"(p)B A (p, m) + 1 > . (7.34) 

For Pauli-Villars's type regularization the function B A (p,m) is decreasing so 
that this condition is implied by the condition at zero momentum 



detM<0 <= 2V"(p)B A (0,m) + 1 > 0. 



For m small we use equation (|4.36 ) and at leading order the condition becomes: 

C(d)(d - 2)m d - 4 V"(p) + 1 > . 



This condition is satisfied by a normal critical point since V"(p c ) > 0. For a 



multicritical point, and taking into account equation ( |7.25|) we find: 



{-l) n C n - 1 {d)m< d -^- d V^ n) {p c ) + 1 > 0. (7.35) 

We obtain a result consistent with our previous analysis: For n even it is al- 
ways satisfied. For n odd it is always satisfied above the critical dimension and 
never below. At the upper-critical dimension we find a condition on the value 
of V^(p c ) which we recognize to be identical to condition ( [7.31[) because then 
2/(ra-l) = d-2. 

The ordered phase. Now m 2 = and the determinant A of the complete 
matrix is: 



-A>0 <S> 2V"(p)B A (p,0)p 2 +p 2 + 4V"(p)a 2 > 0. (7.36) 

We recognize a sum of positive quantities, and the condition is always satisfied. 
Therefore in the case where there is a competition with a disordered saddle point 
only the ordered one can be stable. 

7.4 The scalar bound state 

In this section we study the limit of stability in the disordered phase (a = 0). 
This is a problem which only arises when n is odd, the first case being provided 
by the tri critical point. 

The mass-matrix has then a zero eigenvalue which corresponds to the appear- 
ance of a new massless excitation other than <p. Let us denote by M the p, m 2 
2x2 submatrix. Then 



detM = <(=> 2V"(p)B A (0,m) + 1 = 0. 
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In the two-space the corresponding eigenvector has components (|, V"(p)). 

The small mass m region. In the small m limit the equation can be rewritten 
in terms of the constant C(d) defined in ( |4.16| ): 



C(d) (d - 2)m d - i V"{p) + 1 = 0. (7.37) 

Equation ( |7.37|) tells us that V"(p) must be small. We are thus close to a 
multicritical point. Using the result of the stability analysis we obtain 

(-l) n ~ l d ~ 2 , C n - 1 {d)m n{ - d -V- d V { - n \p c ) = 1 . (7.38) 
(n — 2)! 

We immediately notice that this equation has solutions only for n(d — 2) = d, i.e. 
at the critical dimension. The compatibility then fixes the value of V^ n '(p c ). We 
again find the point ( [7.32| ), V^ n '(p c ) = O c . If we take into account the leading 
correction to the small m behavior we find instead: 

W»)(p c R-'-l~(2„-3)^(^) 4 - J . (7.39) 

This means that when a(d) > there exists a small region V (n) (p c ) > O c where 
the vector field is massive with a small mass m and the bound-state massless. 
The value fi c is a fixed point value. 

The scalar field at small mass. We want to extend the analysis to a situation 
where the scalar field has a small but non-vanishing mass M and m is still small. 
The goal is in particular to explore the neighbourhood of the special point ( |7.32j ) . 
Then the vanishing of the determinant of M implies 

1 + 2V"(p)B A (iM, m) = . (7.40) 

Because M and m are small, this equation still implies that p is close to a point 
p c where V"(p) vanishes. Since reality imposes M < 2m, it is easy to verify that 
this equation has also solutions for only the critical dimension. Then 

VW{Pc)f{m/M) = n c , (7.41) 

where we have set: 

f(z)= / 1 dx[l + (x 2 -l)/(4z 2 )] d/2 " 2 , \<z. (7.42) 
In three dimensions it reduces to: 

f{z) = 2 In 



2z- 1 
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f(z) is a decreasing function which diverges for z = % because d < 3. Thus we 
find solutions in the whole region < V^ n '(p c ) < fi c , i.e. when the multicritical 
point is locally stable. 

Let us calculate the propagator near the pole. We find the matrix A 



A G 2 

where we have set 



dB A (p,m) 



dp 2 



-M 2 



-1 



' n °,2 ) , (7-43) 



p 2 + M 2 \G G 



2 (_Qn-2 W (n) ^ 



G = 7 — : m 

(n-2)! 



For m/M fixed the residue goes to zero with m as m d ~ 2 because the derivative 
of B is of the order of m d ~ 6 . Thus the bound-state decouples on the multicritical 
line. 

The scalar massless excitation: general situation. Up to now we have explored 
only the case where both the scalar field and the vector field propagate. Let us 
now relax the latter condition, and examine what happens when m is no longer 
small. The condition M = then reads 

2V"(p s )B A (0,m) + l = 

together with 

1 f A d d k 

An obvious remark is: there exist solutions only for V"(p s ) < 0, and therefore 
the ordinary critical line can never be approached. In terms of the function F(z) 

F(m /A ) = J k2 + m2 = ^-p J k 2 D{k 2 ) + m 2 ( 7 - 45 ) 

and thus 

r k d ~ 1 dk 

F M= N <lmmr;- (7 - 46) 

the equations can be rewritten 

p s = A d ~ 2 F(z), z = 2V'( Ps )A- 2 , 2A d ~ 4 V" (p s )F'(z) = 1 . 

The function F(z) in Pauli-Villars's regularization is a decreasing function. In 
the same way —F'{z) is a positive decreasing function. 

The third equation is the condition for the two curves corresponding to the 
two first ones become tangent. For any value of z we can find potentials and thus 
solutions. Let us call z s such a value and specialize to cubic potentials. Then 

p s = A d - 2 F(z s ) , 

V(p) = V'(p s )(p - Ps ) + \V"{ Ps ){p - p s ) 2 + ±V^(p s )(p - p s ) 3 , (7.47) 

which yields a two parameter family of solutions. For z small we see that for 
d < 4 the potential becomes proportional to (p — p c ) 3 . 
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7.5 Stability and double scaling limit 

In order to discuss in more details the stability issue and the double scaling limit 
we now construct the effective action for the scalar bound state. We consider 
first only the massless case. We only need the action in the IR limit, and in this 
limit we can integrate out the vector field and the second massive eigenmode. 

Integration over the massive modes. As we have already explained in section 
[7.2| we can integrate over one of the fields, the second being fixed, and we need 
only the result at leading order. Therefore we replace in the functional integral 



e z 



J [dpdX] exp -f trln(-a 2 + X) + N J d d x (-V(p) + \pX) 



(7.48) 



one of the fields by the solution of the field equation. It is useful to first discuss 
the effective potential of the massless mode. This requires calculating the action 
only for constant fields. It is then simpler to eliminate A. We assume in this 
section that m is small (the vector propagates) . For A A the A-equation reads 
(d<4) 

p -p c = -C(d)A (d " 2) / 2 . (7.49) 
It follows that the resulting potential W(p), obtained from equation ([7.171) is 

In the sense of the double scaling limit the criticality conditions are 

w(p) = o((p- Ps r). 

It follows 

V^(p s ) = -\C^[d) ^2/^-2)) rn d ~^ 1 < k < n - 1 . 
For the potential V of minimal degree we find 



W(p) 



r(-2/(d-2)) 



The double scaling limit. We recall here that quite generally one verifies that 
a non-trivial double scaling limit may exist only if the resulting field theory of 
the massless mode is super-renormalizable, i.e. below its upper-critical dimension 
d = In I (n— 2), because perturbation theory has to be IR divergent. Equivalently, 
to eliminate from the critical theory, one has to rescale 

p-p s ocA" 2 V, x^xN {n - 2)e with 1/0 = 2n - d(n -2), 
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where 9 has to be positive. 

We now specialize to dimension three, since d < 3 has already been examined, 
and the expressions above are valid only for d < 4. The normal critical point 
(n = 3), which leads to a </? 3 field theory, and can be obtained for a quadratic 
potential V(p) (the (4> 2 ) 2 ) has been discussed elsewhere. We thus concentrate 
on the next critical point n = 4 where the minimal potential has degree three. 

The d — 3 tricritical point. The potential W(p) then becomes 

W(p)=V(p) + ^(p c -pf. (7.51) 

If the potential V{p) has degree larger than three, we obtain after a local expan- 
sion and a rescaling of fields, 



(^-)(A-m 2 )cxv9/iV, x^Nx, (7.52) 
oZtv p c 



a simple super-renormalizable <£> 4 (x) field theory. If we insist instead that the 
initial theory should be renormalizable, then we remain with only one candi- 
date, the renormalizable (</> 2 ) 3 field theory, also relevant for the tricritical phase 
transition with O(N) symmetry breaking. Inspection of the potential W{p) im- 
mediately shows a remarkable feature: Because the term added to V{p) is itself 
a polynomial of degree three, the critical conditions lead to a potential W(ip) 
which vanishes identically. This result reflects the property that the two saddle 
point equation (dS/dp = , dS/dX = in equations ( [7.24| )) are proportional 
and thus have a continuous one-parameter family of solutions. This results in a 
flat effective potential for (f(x). The effective action for (p depends only on the 
derivatives of <p, like in the 0(2) non-linear a model. 

We conclude that no non-trivial double scaling limit can be obtained in this 
way. In three dimensions with a (4> 2 ) 3 interaction we can generate at most a 
normal critical point n = 3, but then a simple (4> 2 ) 2 field theory suffices. 

The ambiguity of the sign of a(d) discussed in section [Fl| has an interest- 
ing appearance in d = 3 in the small m 2 region. If one keeps the extra term 
proportional to a(d) in equation (|7.50| ) we have 



W(p) = V(p) + 8 -^( Pc - pf + ^4n 2 (p c - pY 



Using now equation ( 7.49 ) and, as mentioned in section [7.4| , the fact that in the 



small m 2 region the potential is proportional to (p — p c ) 3 we can solve for m 2 . 
Since m 2 > the appearance of a phase with small mass depends on the sign of 
a(d). Clearly this shows a non-commutativity of the limits of m 2 / A 2 — > and 
iV — > oo. The small m 2 phase can be reached by a special tuning and cannot be 
reached with an improper sign of a(d). Calculated in this way, m 2 can be made 
proportional to the deviation of the coefficient of p 3 in V(p) from its critical 
value 16-7T 2 . 
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7.6 Conclusions 

This is a study of several subtleties in the phase structure of O(N) vector 
models around multicritical points of odd and even orders. One of the main 
topics is the understanding of the multicritical behavior of these models at their 
critical dimensions and the effective field theory of the 0(iV)-singlet bound state 
obtained in the N — > oo, g — > g c correlated limit. It is pointed out that the 
integration over massive O(N) singlet modes is essential in order to extract the 
correct effective field theory of the small mass scalar excitation. After performing 
this integration, it has been established here that the double scaling limit of (4> 2 ) K 
vector model in its critical dimension d = 2K/ (K — 1) results in a theory of a free 
massless O(N) singlet bound state. This fact is a consequence of the existence of 
fiat directions at the scale invariant multicritical point in the effective action. In 
contrast to the case d < 2K/ (K — 1) where IR singularities compensate powers of 
1/N in the double scaling limit, at d = 2K/ (K— 1) there is no such compensation 
and only a noninteracting effective field theory of the massless bound state is left. 



Another interesting issue in this study is the ambiguity of the sign of a(d). 
The coefficient of m 2 A d_4 denoted by a(d) in the expansion of the gap equation 
in equations ( f7.24cj) and (|7.25|) seems to have a surprisingly important role in 



the approach to the continuum limit (A 2 ^> m 2 ). The existence of an IR fixed 
point at g ~ 0(iV -1 ), as seen in the j3 function for the unrenormalized coupling 
constant (section |4.3|) , depends on the sign of a(d). Moreover, the existence of 
a phase with a small mass m for the O(N) vector quanta and a massless O(N) 
scalar depends also on the sign of a(d). It may very well be that the importance 
of the sign of a(d) is a mere reflection of the limited coupling constant space 
used to described the model. This is left here as an open question that deserves 
a detailed renormalization group or lattice simulation study in the future. 
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